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ABSTRACTS 


The dynamics of a rigid body having two degrees of freedom is 
presented. The complete solutions for nodal points and 
coupled frequencies are obtained. The general solution for the 
nodal positions is presented in graphic form, which provides not 
only a physical picture of any system but also a rapid means of 
evaluating the position of the nodes. Charts have been prepared 
from which the uncoupled frequencies, as required in flutter 
analysis, may be rapidly evaluated from data observed during 
vibration tests. Neglect of this item may produce an unconserva- 
tive error of 20 per cent in the calculated flutter speed. 

A theoretic analysis is made of the factors affecting the power 
and unbalance required to conduct vibration tests. It is shown 
that for a given amplitude the power absorbed by a wing is 
roughly independent of airplane size; however, for any given air- 
plane it is from two to sixteen times greater at torsional than 
bending resonance. At torsional resonance, with one-tenth-in, 
total motion of the trailing edge, a conventional wing absorbs 
about 0.0134 hp. Thus, the theoretic power required is negligible, 
and the choice of a suitable motor depends upon the friction 
horsepower absorbed by the test equipment. 

Numerical examples are presented to illustrate the use of all 
principal equations and charts. 

A discussion is given of the new Curtiss-Wright Universal Vi- 
brator, outlining its major points of interest. 


INTRODUCTION 


7 PRINCIPAL PURPOSE of this paper is to present 
some of the methods that have been developed by 
the St. Louis Plant of the Curtiss-Wright Corporation 
during vibration and flutter analyses made in the past 
two years. The fundamental theory covering the mo- 
tion of a body with two coupled degrees of freedom 
has been derived by several writers.) *% 56 How- 
ever, with regard to coupled frequencies and nodal 
points, experience has shown that the results obtained 
in these references are in a form that cannot easily be 
applied to practical engineering problems. The dynam- 
ics of a rigid body (Fig. 1) with two degrees of freedom 
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Fic. 1. System with two degrees of freedom. Note: X and 
X, are scalar quantities and are always considered positive. All 
other dimensions are shown in positive sense with respect to the 
elastic axis. 


has been derived herewith in a form which, it is be- 
lieved, lends itself most readily to the solution of cer- 
tain vibration problems that often confront the prac- 
ticing engineer. It has been assumed throughout that 
the conventional flutter parameters—namely, elastic 
axis and c.g. location—as well as radius of gyration, 
have been calculated and, in addition, that the 
coupled frequencies have been determined experimen- 
tally. 

The analyses that follow apply rigorously to a rigid 
body supported on springs and having two degrees of 
freedom. As discussed later, it has been assumed 
throughout that the results obtained from such a system 
can be applied to actual three-dimensional structures in 
which the elasticity is distributed throughout the struc- 
ture itself. Experience, supported by Theodorsen‘ 
and several other investigators,* ® has justified this 
assumption. In many cases of airplane vibrations the 
results can be applied quantitatively with a high degree 
of accuracy, and in almost every case they can be 
applied qualitatively. 

Solutions to equations have been presented graphi- 
cally to facilitate practical engineering work, and nu- 
merical examples have been worked out to illustrate the 
applications of the principal charts and equa- 
tions. 
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has been derived herewith in a form which, it is be- 
lieved, lends itself most readily to the solution of cer- 
tain vibration problems that often confront the prac- 
ticing engineer. It has been assumed throughout that 
the conventional flutter parameters—namely, elastic 
axis and c.g. location—as well as radius of gyration, 
have been calculated and, in addition, that the 
coupled frequencies have been determined experimen- 
tally. 

The analyses that follow apply rigorously to a rigid 
body supported on springs and having two degrees of 
freedom. As discussed later, it has been assumed 
throughout that the results obtained from such a system 
can be applied to actual three-dimensional structures in 
which the elasticity is distributed throughout the struc- 
ture itself. Experience, supported by Theodorsen* 
and several other investigators,®*® has justified this 
assumption. In many cases of airplane vibrations the 
results can be applied quantitatively with a high degree 
of accuracy, and in almost every case they can be 
applied qualitatively. 
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SYMBOLS 


(All symbols are in fundamental ft.lb. sec. units) 





*a = elastic axis parameter 
*b = half chord of wing 
B = wing span ° 
" = bending stiffness coefficient per unit span of two- 
dimensional wing 
"C. = torsional stiffness coefficient per unit span of two- 
dimensional wing 
F = force applied by vibrator 
g = acceleration of gravity 
* 205 Zn = damping coefficients 
‘ = deflection of elastic axis 
hy = deflection of point of application of F 
hy = deflection of any reference point 
L = *bxg = distance between c.g. and elastic axis 
L, = distance between rotors of a twin rotor vibrator 
a = power 
r = *bre = radius of gyration about elastic axis 
W = weight of equivalent two-dimensional body 
Wore = unbalance of vibrator 
Xr = distance from node toreference point (always posi- 
tive) 
Ag = distance from node to vibrator (always positive) 
a = angular deflection (rad.) 
ri) = phase angle (rad. 
My Me, Mn = see definition following Eq. (19 
"x = mass ratio 
*y = frequency of applied force (rad. per sec.) 
2 = frequency of applied force at resonance (rad. per sec.) 
Wy, Ws = greater and smaller uncoupled frequencies of sys- 
tem with two degrees of freedom (rad. per sec.) 
"Wiss = uncoupled torsional frequency about elastic axis 
(rad. per sec.) 
*wn = uncoupled bending frequency (rad. per sec.) 
29, Qs = greater and smaller coupled frequencies of system 
with two degrees of freedom (rad. per sec.) 
unbalance of smaller vic ; : 
r = oi for a twin rotor vibra- 
( unbalance of larger rotor 


tor (positive when rotors are in phase; negative 
when 180 deg. out of phase 

bL, = distance from point midway between rotors to 
single equivalent rotor of a twin-rotor vibrator 


= distance from elastic axis to applied force (positive 


gr - 
when elastic axis lies between c.g. and force) 

cr = distance from elastic axis to node (positive when 
elastic axis lies between c.g. and node) 

v1, fs = positive and negative solutions, respectively, for ¢ 


at resonance 


THEORY 

Reference is made to the system illustrated in Fig. | 
which represents a rigid body having two degrees oi 
freedom. It is supported at the elastic axis by means of 
a spring of stiffness C,, while restraint about the elastic 
axis is supplied by a coil spring of torsional stiffness C,. 
An oscillating force, F = F) cos wt, is applied at a dis 
tance gr from the elastic axis. The equations of motion 
may be determined by equating to zero the sum of all 
forces and moments. Whence 


2F=0 

— (d*h/dt®? + Ld*a/dt?)(W/g) —hC, + F=0 (1) 
=i =0 (about elastic axis) 

— (Ld*h/dt? + r°d?a/dt?)(W/g) — aC, — Fir =0 (2) 


Assuming sinusoidal motion and making the following 
substitutions 


lw wi 
h = Ie™ a= ae 
d*h/dt? = = whe d*a ‘dt? = wae 
F = Foe 


Ci/(W/g) = wr? C./177(W/g) = we? 


Eqs. (1) and (2) become 
ho(w? — w,”) + ragww?L/r + Fo/(W/g) = 0 (3) 
how*L/r + rao(w? — we?) — Fot/(W/g) = 0 (4) 


It should be noted that Eqs. (3) and (4) define com- 
pletely the steady state motion of the body in response 
to a sinusoidal force of amplitude Fo and having any 
frequency w. 

With reference to Fig. 1 it can be seen that the point 
about which the resultant response occurs at any fre- 
quency w is given bv 


“" 1 
¢ = No/ ra 


At resonance, the nodal points and frequencies are 
given by the particular solution of Eqs. (3) and (4) 
(i.e., for Fy = 0). 

Therefore, when w is equal to the resonant frequency 
Q and ¢ is substituted for ho/rao, Eqs. (3) and (4) become 











: es , i ‘ *(Q2 — ,.2) "| (5) 
Subscript 0 indicates amplitude of oscillating function. g(Q w,*) + Q L/ r=0 \O) 
* These symbols are the same as those used by Theodorsen fQ2L/r + (Q? — w,”) = 0 (6) 
whence 
” (wWe,/‘wn,)? — 1] += V [(con/wn)? — 1]*? + 4(we/w,)?(L/r)? J 
f=— —— G 
2L/r 
i. is 2 ~/ 2 2)2 ae r 2 / 2) 

o2 (We" + wa”) = (Wa? + wp)? — 4eQ7w,2(1 — L?/r?) (8) 

LY a => — c 





2(1 — L?/r*) 


By inspection of Eq. (7) it can be seen that the sign before the radical always determines the sign of {, the positive 
value being defined as ¢; and the negative by f. Therefore there are always two nodes—one on each side of the 
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elastic axis. From Fig. 1 it can be seen that the nodal points lie outside of the range between the c.g. and elastic 
axis for all values of {, and for values of {2 numerically greater than L/r. For any given value of L/r, {2 decreases 
numerically as w/w, increases. Thus, the smallest possible numerical value of {; is obtained when w,/w, ap- 


proaches infinity. Eq. (7) may be written 


f= (wWe./w,)? — ; tes AP 4+ 4 (a/eon)"(L/r)*_ / 


2L/r [(wa/on)®? — 1]? | 








As (wa/w,) becomes very large (we/w,)? — 1 approaches (w,/w,)?. Expanding by the binomial theorem and re- 
taining only infinitesimals of the first order, the expression for {2 may be written 


Be b-f)) 


lim £2 (coee/en) —-> oS L/r (9) 


Finally 


Therefore the negative root can never be numerically 
less than L/r. Thus it can be seen that the nodal 
points can never lie between the c.g. and the elastic 
axis. 


Again from Eq. (7) 
fife = — (wa/wn)? (10) 


In order to identify the corresponding values of ¢ and w, 
Eq. (5) may be written 


(Q/wy)? = ¢/(¢ + L/r) (11) 





By inspection of Eqs. (9) and (11) it can be seen that 
¢/(¢ + L/r) is always positive. Further, from Eq. (11), 
(Q/w,)* is greater than unity for all negative values of ¢ 
and less than unity for all positive values. Thus the 
positive root of Eq. (7) is always associated with the 
smaller value of 2 and vice versa. Alternatively, it 
may be said that the larger coupled frequency Q, always 
occurs about a node that lies on the same side of the 





elastic axis as the c.g. Conversely, the c.g. and the 

node for the lower frequency , lie on opposite sides of 

the elastic axis. Eq. (7) has been presented graphically Fic. 2. Location of nodal points at resonance. (For reference 
ae ot 9 9 see Eq. (7).) Note: (a) and (b) are identical except for the,scale 
in Figs. 2a and 2b. of f. 


The foregoing discussion has applied only to the 
nodes that occur during free vibrations at resonance. 
Eqs. (3) and (4) may be solved to find 4p and rag. The sponse to an applied force of any frequency w is found to 
quotient of these solutions /p/ray is equal to ¢, from be . 


which the general solution for the nodal point in re- 








f= { lw? — we?/(EL/r + 1)]/[o? — fon?/(L/r + SS L[E(L/r) + 1/[(L/r) + 8} (12) 
From Eq. (12) it can be seen that during forced vibrations ¢ becomes independent of the applied frequency w if 
wa®/[E(L/r) + 1] = kon?/[(L/r) + (13) 
Eq. (13) may be solved for ¢, hence 
¢ — [ea/en)? = 1] + VUwa/en)? = 1]? + 4(a/ eon) (L/ 0) (14) 
2L/r 


Comparison of Eqs. (7) and (14), which give identical expressions for § and {, shows that when the vibrator is 
placed at one node the point about which resultant rotation occurs is independent of the applied frequency and 


must be coincident with the other node. 
Inspection of Eq. (8) shows that, as far as the coupled frequency 2 is concerned, wa and w, are interchangeable. 
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For this reason the uncoupled frequencies should be replaced by w, and w;, where the subscripts indicate “‘large’’ 
and ‘‘small’’ instead of “‘torsion’’ and “‘bending.’’ Making these substitutions in Eq. (8) it is found that 














2,2 + 2,2 = (wi? + w,?)/(1 — L*/r*) (15) 
0,72,? = wy?w,?/(1 — L?/r?) (16) 
From Eqs. (15) and (16) it is found that 
lave 1 — L*/r*)(1 + 2,2/0,2) + V (1 — L?/r?)2(1 + 2,2/2,2)? — 4(1 — L?/r2)0,2/2,2 a 
(wy/ Q,)? = > mae (1 7) 
fo \2 . A — L*/r)(1 + 2,2/2,2) — VL = L*/P)1 + 2,2/2,%)? — 41 — L*/r)2,2/2,? a 
(w5/ 2,)° = 2 (18) 


These equations have been presented in graphic form in Fig. 3, from which the uncoupled frequencies may be 
found corresponding to any known values of L/r and the coupled frequency ratio 2,/Q,. It will be noted that the 
information given by Eqs. (17) and (18) is insufficient to determine which of the coupled frequencies is torsion 
about the elastic axis and which is bending. This is due to the fact that in Eq. (8) wa and w, are interchangeable. 
However, this is not true of Eq. (7), which can be used to establish the identity of the uncoupled frequencies. 
It can be seen that Eq. (7) would have been sufficient to find w./w, directly if the values of ¢ and L/r had been 
known. However, in vibration testing it is difficult to determine the nodal points with any degree of accuracy, 
whereas the measurement of frequency is quite straightforward. In this way it is most expedient to use Eq. (8) 
(or Fig. 3) to find the uncoupled frequencies, while a rough check between the observed nodes and those estimated 
by means of Eq. (7) (or Figs. 2a and 2b) will usually remove all ambiguity as to the identity of wz and w,,. 

The previous discussion refers strictly to a body vibrating in a vacuum. The equations of motion for a flat 
plate vibrating in the atmosphere have been derived by Theodorsen.* At zero air speed and assuming the contro] 
surface to be rigid, the solution of Eqs. (A) and (C) on page 10 of reference 3 reduces to the form 


o(y — J*2 b — salad” a ~ + 
ra®{] (1 + «)[1 + «(/e + a°)/ret (I+ 4%) © [L+«(/s+a*)/re'] 


con® se =0 (19) 
(+ thi + «(Ys + ay 

















Introducing the following parameters, which all become unity for x = 0 (i.e., in a vacuum), 








ut = (1 — ax/xa)*/{(1 + «)[1 + x("/s + a*)/re']} 
pa? = 1/[1 + w(1/s + @*)/ra*] 
~ = 1/(1 k) 
Eq. (19) takes the form 
2 = Muaita)® + (uneon)*] + V [(uatra)® + (uneon)*]* = 4(uatra) *(unton)*[1 = (usta/ra)I (20) 


2[1 — (uxa/ra)*| 


It can be seen that Eqs. (8) and (20) express exactly the same functional relationship if (ux./Ta), (Mata) and (upp) 
be substituted for (L/r), wa and w,, respectively. The parameters y. and yu, cannot be arbitrarily associated with 
the greater and smaller frequencies as was done with w, and w, in Eqs. (17) and (18). However, the air mass cor- 
rection can be applied only to bodies having the form of flat plates or airfoils; further, it may be said that for all 
conventional airfoil surfaces w, is greater than w,. Thus the solution for the uncoupled frequencies, as given by 
Eq. (20), has been determined in the form of Eqs. (17) and (18), except that it has been assumed that paw, is al- 
ways greater than y,w,. These results have been plotted in Fig. 4. 


APPLICATION TO THREE-DIMENSIONAL STRUCTURES 


All modes of vibration exhibited by an airplane are the whole airplane, this assumption is well justified. 
functions of the overall elastic and inertia character- Even so, a complete analysis of each three-dimensional 
istics of the structure. For the sake of simplicity, it is component is usually complicated and may be replaced 
common practice to consider the major components by a two-dimensional system. The choice of the 
separately—for example, a wing, an engine mount, the equivalent system depends entirely upon the ingenuity 
tail unit, etc. Provided that the mass and stiffness of of the designer; however, several commonly used 
each component is relatively smal] compared to that of examples will be discussed. 
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6 1.0 


O 2 4 6 ; 
= 1/1, — : 


Fic. 3. Uncoupled frequencies ina vacuum. (For reference see 
Eqs. (17) and (18).) 


Wings 

The equivalent weight may be taken as one-sixth of 
the entire weight of both wings and center section or, 
alternatively, as the weight of the outer half of one 
panel. The latter choice is better but frequently can- 
not be easily obtained. 

The equivalent radius of gyration, c.g. and elastic 
axis locations may be taken to be identical with those of 
a section chosen three-fourths of the semispan from the 
centerline of the airplane. 

The amplitudes of the equivalent system will also cor- 
respond to those of the actual wing at this reference sec- 
tion. Roughly speaking it may be said that the tip 
amplitude is twice the amplitude of the equivalent sys- 
tem—.e., of the reference section. 


Tail Units 


With respect to fuselage torsion and side bending, the 
entire tail unit may be considered as a rigid body, its 
weight, c.g. location and radius of gyration being used 
directly in the equivalent system. The actual elastic 
axis of the fuselage determined at the rudder hinge sta- 
tion supplies the remaining parameter. 


Engine Installations 


The engine itself may be considered as the rigid body. 
Provided that one of the principal elastic axes of the 
mounting structure is nearly parallel to the thrust 
line, it may be used in the equivalent system. If, 


(HnWn /Ms) 


——(PaWe/ne) 





+ 6 8 1.0 
=N,/N, es 


a 


Fic. 4. Uncoupled frequencies in air. (For reference see Eq. 
9 


«“ 


however, it is not approximately parallel, this problem 
must be given special consideration. 


APPLICATION TO FLUTTER ANALYSIS 


In the fundamental flutter theory, the following pa- 
rameters, as defined by Theodorsen,*‘ are of primary im- 
portance: w, = uncoupled circular bending frequency 
(rad. per sec.) and w, = uncoupled circular torsional 
frequency about the elastic axis (rad. per sec.). These 
parameters may be calculated, but more often they are 
determined by means of vibration tests on the finished 
airplane. Frequently, the observed coupled frequencies 
are substituted directly for w, and w,. In general, the 
elastic axis and c.g. of a wing are not coincident. For 
this reason the observed coupled frequencies are not 
identical with the desired uncoupled frequencies, as 
shown in Eq. (8). In many cases the error involved is 
not great, but, as the coupling factor L/r increases or 
the uncoupled frequency ratio approaches unity, the 
discrepancy becomes of increasing importance. For 
the purpose of analyzing vibration tests Eq. (8) has 
been written in the form of Eqs. (17) and (18). These 
have been plotted in Fig. 3, by means of which the un- 
coupled frequencies may be immediately determined 
from the observed values of 2, and Q, and the calcu- 
lated value of L/r. It must be noted that L/r is nu- 
merically equal to x./fa, as defined by Theodorsen in 
his flutter theory.‘ The effect of the surrounding atmos- 
phere may be taken into account concurrently with the 
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correction for coupling by means of Eq. (20). The solu- 
tion has been plotted in Fig. 4. 

It must be noted that all parameters required for 
these transformations must be found in any case during 
a routine flutter analysis. Consequently, the addi- 
tional work entailed by these corrections is quite 
negligible. 


APPLICATION TO VIBRATION TESTING 


There are two principal purposes of vibration testing: 
(1) to determine the resonant frequencies of the test 
specimen in its fundamental modes and (2) to deter- 
mine the position of the nodal points corresponding to 
the fundamental modes. 

Difficulty is encountered under the following condi- 
tions: (a) When the coupled frequency ratio is near 
unity it becomes difficult to distinguish between the 
two modes and, in particular, the nodal points become 
indistinct due to the superposition of one mode on the 
other. This condition is not uncommon with respect to 
fuselage bending and torsion. (b) When extraneous 
responses (e.g., skin panel buckling or harmonics) oc- 
cur, the identification of the fundamental modes some- 
times presents quite a problem. 

In case (a) the problem of separating the two funda- 
mental modes may be greatly simplified by placing the 
vibrator at one node. When this has been done, only 
vibrations about the other node can be excited, as 
shown by Eqs. (7) and (14); in this way each mode can 
be investigated separately without interference from 
the other. The Curtiss-Wright Universal Vibrator, as 
will be described later, is particularly well adapted to 
making use of this phenomenon. 

In case (b) the extraneous responses may have widely 
differing characteristics. At times there is no doubt 
about their identification, while at others they become 
most confusing. In either event the observer can rely 
on his experience and judgment to make the proper 
choice. <A careful comparison between the observed 
motions and the anticipated characteristics as indicated 
by Figs. 2a and 2b will often yield a definite answer to 
an otherwise ambiguous problem. Specific rules cannot 
be given because of the widely differing types of local 
responses that may be observed. 

Even though a certain node may be well defined dur- 
ing a test, it is often extremely useful to have avail- 
able some ready theoretic means of substantiating its 
position. That this is true can best be illustrated by 
means of an actual example. During the vibration 
tests of Airplane B (see ‘‘Numerical Example’’) the 
fuselage torsional node was distinctly observed to lie 
approximately 1 ft. below the tail cone. Because of the 
fact that the dynamic balance coefficient of the rudder 
nad to be calculated with respect to this node, it was 
important that it should be observed correctly. A 
hearty discussion arose among the various observers, 
most of whom contended that it could not possibly lie 
outside of the fuselage and that probably the whole 


airplane was moving and therefore the node appeared 
in a fictitious position. The numerical result obtained 
by applying the methods of this paper to Airplane B 
shows clearly that not only was it possible for the node 
to lie outside the fuselage but also that the observed 
node was entirely correct. Had the information of 
Figs. 2 and 3 been available at the time of the test, the 
discussion might have reached an earlier and more 
rational conclusion. 


POWER AND UNBALANCE REQUIRED FOR VIBRATION 
TESTS 


The choice and design of vibrators depend largely 
upon the three following items of fundamental impor- 
tance: 

(1) The power required to produce a given ampli- 
tude, the choice of the latter being a function of the 
available means by which it may be observed. When 
sensitive seismic-type pickups are used, an amplitude 
of only a few hundredths of an inch produces a measur 
able oscillation, while for direct visual observation about 
one-tenth of an inch total motion is necessary. 

(2) The unbalance needed to absorb the required 
horsepower and, conversely, to produce the desired 
amplitude. 

(3) The effect of airplane size upon the required 


power and unbalance. 

It has been shown by Theodorsen‘ that the internal 
damping of a wing is in phase with the velocity during 
vibration, while its magnitude is proportional to the 
amplitude of the elastic restoring force. Neglecting 
the effect of air damping, which can be shown to be 
quite negligible at zero air speed, the following may be 
concluded (when h = fy sin wt and a = ay sin wt): 


Work done per cycle by internal damping 


‘x /2 
oa if 2nCrho cos wt dh + LaCaa cos wt da (21) 
0 


‘ae /2 
= 4(g,Crho? + faCucn') cos? wf d(wt) 
0 


= (gnCyho® + gaCaco”)e (22) 
But 
Cy = wn? W/¢ 
Ca = r'w’W/zg 
From Fig. 1 it can be seen that 
No = hy (&/X) 
a = (hw/X,) 


Work done per cycle by damping 
Time of one cycle 





Power = 


99 


Therefore, making the above substitutions in Eq. (22) 


P = 1/o(W/g)hy?w* [gn (wn/w)?2(¢/X)? + 
ga(wa/w)?(1/X)*] (23) 








Sir 
lie 
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The power absorbed by the vibrator will now be 


considered. Assuming sinusoidal motion 
h, = hy sin wt 
F = Fy sin (wt + ¢) = (W,r,/g)w* sin (wt + ¢) 
Power absorbed by vibrator 
= JS F dh,/(2r/w) 
= 1/,(W,r,/g)w*h» sin 
At resonance ¢ = 2/2 
P = '/2(Wor,/g)w*lo (24) 
From Fig. 1 
ho = Ayw(X,/X) 


Substituting in Eq. (24) 


Pr == 1/5(W or, /g)w*hye(X ,/X ) (25) 
Eliminating P from Eqs. (23) and (25), rearranging 
and introducing the wing span B 
Wre/Whw = (1/Bu)*(X/X,) [ga(Bon)*(¢/X)? + 
a(Bwa)?(1/X)*] (26) 


Comparison of Power Required for Bending and Torsion 


Experience has shown that for conventional wings 
the torsional frequency is usually several times as great 
at that for bending. Thus, at bending resonance it 
can be seen from Figs. 2a and 2b that the nodal point 
lies several chord lengths ahead of the wing, while for 
torsion it lies near the mid chord point. The trailing 
edge will be assumed to be the observation point at 
which the amplitude fy is measured. Then for all 
conventional wings the approximate values given in 
Table 1 may be used. Substituting the average values 
of Table 1 into Eq. (23), and assuming that ge = gp, 
for bending (when w ~ aw») 


P = 1/o(W/g)hw*o,*g, (0.81 ao 0.05) 
Therefore 


P/hy? ~ */2(W/g)g,w,* (for all wings at bending 
resonance) (27) 
Similarly, for torsion (when w ~ w,) 


P = '/2(W/g)hwo*gawa® (0.03 + 0.25) 

Therefore 

P/hy? ~ '/2(W/g)g..*/4 (for all wings at torsional 
resonance) (28) 

A comparative study of Eqs. (27) and (28) shows that 

for the same trailing edge amplitude 





P . 
(at torsional resonance) os 1/4(q/w,)® 
P (at bending resonance) 
Since the frequency ratio of conventional wings usually 
lies between 2 and 4, it can be seen that the power re- 
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TABLE 1 


Range of Values of Parameters for Conventional Wings 





Value at Bending Value at Torsional 


Parameters Resonance Resonance 
o/X 0.8 —1.0 —-0.3 ~—-0.5 
1/X 0.05 + 0.10 0.4—> 0.6 
wp/ w 1.0 0.25—~ 0.50 
We, /w 2—4 1.0 


quired at torsional resonance is from two to sixteen 
times greater than that required to produce the same 
trailing edge amplitude in bending. 

From Eqs. (27) and (28) it can be seen that for ‘a 
given amplitude, which is determined by the available 
measuring equipment, and for geometrically similar 
wings 


P « Ww « (Bw)? 


Experience has shown’ that the resonant frequencies 
vary inversely as the span. Therefore, the power re- 
quired to produce a given amplitude is roughly inde- 
pendent of airplane size. 

It will be shown later by means of a numerical ex- 
ample (see Airplane A) that the theoretic power re- 
quired to produce one-tenth of an inch motion at the 
trailing edge is 0.00064 and 0.0134 hp. at bending and 
torsional resonance, respectively. Since the actual 
power required to vibrate the wing is so small, it can 
be seen that the motor for the vibrator need only be 
designed to overcome the friction horsepower in the 
moving system. 

A similar analysis of Eq. (26) shows that the required 
unbalance varies directly with the amplitude and with 
the weight of the surface being tested. Thus, for a 
given amplitude, the required unbalance must increase 
roughly as the cube of the wing span. 

The foregoing analysis can be applied to a three- 
dimensional wing by means of the substitutions outlined 
under the heading ‘‘Application to Three-Dimensional 
Structures.”’ 


NOTE ON STRUCTURAL DAMPING 


At the present time there is little quantitative infor- 
mation about the structural damping present in aircraft 
structures. A limited amount of experimental data 
has been summarized by Wylie.’ The actual values of 
g, and gq depend upon the type of design and the ma- 
terial used. They vary from 0.04 to 0.09 for metal 
monoplanes, and usually g, ~ ga. For purposes of 
numerical examples in this paper, the structural 
damping has been assumed to be 0.05 in all cases. 
It can be seen from Eqs. (23) and (26) that the power 
and unbalance required to produce a given amplitude 
both vary directly with the structural damping co- 
efficients. 

For a particular structure the damping coefficients 
may be evaluated from test results if the amplitude re- 
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sponse to a constant exciting force at different fre- 
quencies is observed. Then 


g = (m2? — m*)/(m2? + m”) 


where ; and m2 are the two frequencies at which the 
observed response is equal to the maximum response 
divided by the square root of 2. The appropriate sub- 
script must be supplied according to whether bending 
or torsional oscillations are observed. It must be noted 
that a constant applied force is specified. When a con- 
ventional unbalanced vibrator is used, all observed 
amplitudes must be divided by the square of the cor- 
responding frequency before finding 1; and m. 


NUMERICAL EXAMPLE 


Airplane A 


Chord = 4.71 ft. 

r = 1.48 ft. 
Xa = 0.36 

Pex = 0.63 
Coupling factor = L/r = xa/ra = 0.57 

a = — 0.40 

k = 0.12 
Weight of entire wing unit = 730 lbs. 
W = 730/6 = 122 Ibs. 
Measured bending frequency = 580 c.p.m. 
Measured torsional frequency = 2,400 c.p.m. 


(1) Determination of Uncoupled Wing Frequencies 














2, = (580) (27/60) = 60.8 rad. per sec. 
Q, = (2,400) (27/60) = 251.2 rad. per sec. 
Coupled frequency ratio = 2,/Q, = 0.242 
| elgg 7% 2 
= | [1 (2x /*a)) rs = 1.026 
VL + «tl + (k/ra*)(*/s + 2*)} 
Ms 
hea = | . : = 0.960 
V1 + (k/ra*)(*/s + a’) 


0.945 


Il 


wy = V1/(1 + «) 
From Fig. 3 (neglecting effect of air mass) 


= 1.016 
= 0.809 


w,/Qs; 

w,/Q, 
Because of the conventional nature of the wing vibra- 
tion, it is quite evident that w, is associated with w, 
and w, with w,, and there is no ambiguity. 


TABLE 2 


Example of Coupled and Uncoupled Frequencies 





-—Uncoupled Frequencies— 


Measured Air mass Air mass 

Frequency, neglected, corrected, 
c.p.m. c.p.m. c.p.m. 
Bending 580 588 626 
Torsion 2,400 1,940 1,995 
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w@, = (1.016)(580)(27r/60) = (588) (27/60) rad. per sec. 
we = (0.809) (2,400) (2x/60) = (1,940)(2r/60) rad. per 
sec. 


From Fig. 4 (including air mass correction) 


MXa/Ta = (1.026)(0.57) = 0.585 

Myo, /Q, = 1.018 

MaWa/Q, = 0.798 

wr, = 1.0182,/u, = [(1.018)(580) /0.945] [27/60] 
= (626)(27/60) rad. per sec. 

We = 0.7980,/ue = [(0.798) (2,400) /0.960] [27/60] 


= (1,995) (22/60) rad. per sec. 


The results of the numerical example have been sum- 
marized in Table 2. Thus it can be seen that the 
effect of coupling is not negligible insofar as it allows 
nearly a 20 per cent error in the torsional frequency 
that should be used in the flutter calculations; this en- 
tails, roughly speaking, a 20 per cent unconservative 
error in the calculations of the bending-torsion flutter 
speed. The effect of air mass is considerably less, 
being in the order of 5 per cent for both, though some- 
what greater for bending than torsion. Further, the 
neglect of the air mass correction leads to a conserva- 
tive answer. 


(2) Power and Unbalance Required to Vibrate Wing 
From Table 2 
(1,995) (27/60) = 208.8 rad. per sec. 
w, = (626)(27/60) = 65.6 rad. per sec. 
We /Wp = 3.18 


Wa 


Therefore from Figs. 2a and 2b 


= +16.0 
= —0.63 


, Sy 
_ 


wn 
i) 


The elastic axis is at 30 per cent of the chord from the 
leading edge. 

With reference to Fig. 1 it can be seen that at bending 
resonance 


Xr = fir + 0.7 (chord) 
16 + [(0.7)(4.71)/1.48] = 18.2 
$:/X = 0.88 


X = 


Similarly, at torsional resonance 
X = — 0.63 + [(0.7)(4.71)/1.48] = 1.60 
fo/X = — 0.39 


Assuming that gz = g, and substituting the appro- 
priate values in Eq. (23), at bending resonance 


° 3. /65 @\2 
ws $0.84 ($8.8) (0.88)? + 


32.2 2 60.8 


(78)( 1 y] 
60.8/ \18.2 


At torsional resonance 


122 


P= hyo*Zn x 











= hy*g, X 0.40 X 10® (29) 


























GROUND 


95] 93 > 2 
x 1.3 (6) (—0.39)? + 


32.2 2 251.2 


9 2\\ 2 2 
(5) (—) | = Ine X 8.3 X 10° (30) 
251.2/ \1.60 


With reference to these numerical operations it can be 
seen that the power required for bending depends almost 
entirely upon the first term in parentheses, while that 
for torsion is controlled by the second term. Eqs. 
(29) and (30) show that the power required is 20 times 
greater for torsion than for bending if the trailing edge 
amplitude is to be the same in both cases. 

In order to estimate the required unbalance, it shall 
be assumed that the vibrator is located on the trailing 
edge, in which case X,/X = 1. Making the appro- 
priate substitutions in Eq. (26) and canceling out the 
span B (since the actual frequencies are known), at 
bending resonance 





122 








P= hy" Ze x 





Wr,/Whw = gn X 0.937 (31) 


At torsional resonance 


W.r,/Whw = ga X 0.280 (32) 


The power and unbalance required to produce one- 
tenth of an inch total motion at the trailing edge will 
be determined. 


0.05/12 = 0.0042 ft. 


hyo i 


Substituting in Eqs. (29), (30), (81), and (32) and as- 
suming that ga = g, = 0.05 


(0.0042)? K 0.05 X 0.40 X 108 





Horsepower = —- 
550 
= 0.00064 hp. (at bending resonance) 

Horsepower = 0.0134 hp. (at torsional resonance) 
Unbalance = W,r, = 122 X 0.0042 X 0.05 X 

0.937 X 12 = 0.288 Ib.in. (at 

bending resonance) 
Unbalance = 0.086 lb.in. (at torsional resonance) 


The numerical values obtained for the horsepower re- 
quired are roughly representative of all airplanes simi- 
lar in size and general design to the one chosen. This 
is not true, however, with regard to the unbalance re- 
quired; in this case the location of the vibrator has 





Fic. 5. Tail unit of airplane B. 
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such a large effect that each particular location must be 
investigated separately. 


Airplane B (Analysis of fuselage vibrations) 


Vibrations involving torsion and side bending of the 
fuselage shall be considered. The tail end of this air- 
plane is shown diagrammatically in Fig. 5. 


L = 35 in. (estimated) 
r = 83 in. (estimated) 
Coupling factor L/r = 0.42 
Side bending frequency (measured) = 700 c.p.m. 
Torsional frequency (measured) = 350 c.p.m. 


From Fig. 3 


w,/Q, = 1.05 
w,/2, = 0.88 
1.05 350 
w,/o, = —— X ~ = 0.60 
0.88 700 


It is now known that the uncoupled frequency ratio is 
I | : 
0.60, but it is uncertain where we or w, is the larger. ° 
With reference to Figs. 2a and 2b, assuming 


We/a, = 1/0.60 = 1.7 
ie = +5.5 
Se as —0.57 


The nodal points are given by {r—i.e., 460 in. below 
the elastic axis and 47 in. above the elastic axis. These 
do not agree with observed results. Therefore it shall 


be assumed that 
We/w, = 0.60 
from Figs. 2a and 2b 


o = 0.20 
Se == —1.8 


Hence, the nodal points are 17 in. below elastic axis 
and 150 in. above elastic axis. These agree with the 
observed results as discussed after Eq. (18) and also 
at the end of the discussion on “Application to Vibra- 
tion Testing.” 


CURTISS-WRIGHT UNIVERSAL VIBRATOR 


The vibrator is illustrated diagrammatically in Fig. 
6, while photographs are shown in Figs. 7 and 8. 
Fundamentally it consists of two unbalanced rotors 
mounted on parallel shafts 20 in. apart. They are 
driven in the same direction by means of bevel gears 
mounted on a stiff aluminum alloy torque tube. In 
this way the natural frequency of the system is above 
the highest speed for which the vibrator has been de- 
signed to operate (viz., 2,500 c.p.m.). The rotors may 
be connected so that they are in phase (i.e., so that cen- 
trifugal forces act in same direction) or 180 deg. out of 
phase as shown in Fig. 6. In either case the combination 
can be replaced by a single equivalent rotor at a point 
R, as shown in Fig. 6b. When the rotors are in phase, 









490 JOURNAL OF THE AERONAUTICAL -SCIENCES—NOVEMBER, 1942 





eo (a) 
PURE COUPLE ‘e 


oe.) a 


Qe Ay 














benign 
. UNBAL (© SINGLE 
~ UNBAL. @ JE QIV T_ 
(ROTOR 
-t I-dr 
8 =a Fx 
(b) 


Fic. 6. Curtiss-Wright Universal Vibrator (diagrammatic). 


R lies between them and the equivalent unbalance is 
equal to the sum of the actual unbalances; when the 
rotors are 180 deg. out of phase, R lies to one side as 
shown in Fig. 6b and the equivalent unbalance is equal 
to the difference between the unbalances of the actual 
rotors. The location of R with respect to the center- 
line of the system is defined by 6L,, where L, is the dis- 
tance between rotors. The parameter 6 is given by the 
equation 


6 = '/(1 — A)/(1 + A) (33) 
where X is the ratio of the smaller to the larger rotor 
unbalance. When the rotors are in phase, \ is positive; 
when 180 deg. out of phase, \ is negative. From Eq. 
(33) it can be seen that the single equivalent rotor 
may be made to act at any point along the line joining 
the two actual rotors. In the limit, when the rotors 
are equal and A = —1, the equivalent rotor becomes 
infinitely small and is located at infinity. In this way 
a pure couple is applied by the vibrator, which is par- 
ticularly effective for the purpose of exciting torsional 
oscillations. 





[eae 


ae 





Fic. 7. Curtiss-Wright Universal Vibrator. 





Fic. 8. Detail of rotor. 


In order to make use of the phenomenon illustrated 
by Eq. (33), the unbalanced weight on each rotor is 
built up by means of eight separate sectors cut from 
1/s-in. steel sheet. In this way the unbalance of each 
rotor may be varied from 1 to 7 lb.in. As explained 
previously this feature can be used to advantage in iso 
lating obscure modes. 


REFERENCES 


* Den Hartog, J. P., Mechanical Vibrations, Second Edition, 
p. 102; McGraw-Hill Book Co., Inc., New York, 1940. 

2 Timoshenko, S., Vibration Problems in Engineering, Second 
Edition, p. 232; D. Van Nostrand Co., Inc., New York, 1937. 

3 Theodorsen, Theodore, General Theory of Aerodynamic In- 
stability and the Mechanism of Flutter, N.A.C.A, T.R. No. 496, 
1935 

* Theodorsen, Theodore, and Garrick, I. E., Mechanism of 
Flutter. A Theoretical and Experimental Investigation of the Flutter 
Problem, N.A.C.A. T.R. No. 685, 1940. 

5 Lombard, A. E., Flutter in Aircraft, Thesis in partial fulfill- 
ment of the Degree of Doctor of Philosophy, California Institute 
of Technology, Pasadena, Calif., 1939 (unpublished). 

® Kassner, R., and Fingado, H., The Two-Dimensional Problem 
of Wing Vibration, Journal of the Royal Aeronautical Society, 
Vol. 41, pp. 921-944, October, 1937; translated from Luftfahrt- 
forschung, Vol. 13, No. 11, pp. 374-387, November 20, 1936. 

7 Wylie, Jean, Perpendicular Axes Control Surface Binary 
Flutter, p. 184; Civil Aeronautics Administration. Certificate 
and Inspection Release No. 46, Aircraft Airworthiness Section 
Report No. 23, April, 1941. 








fe - ses te ome 


ol 


acl 


sior 














Determination of Airplane Critical Altitude 
by Flight Tests 


L. C. MILLER* 


Brewster Aeronautical Corporation 


SUMMARY 


The determination of airplane critical altitude has become of 
considerable importance during the last four or five years. The 
method presented in this paper has been developed around a 
practical problem on a specific airplane. By use of this method 
it has been possible to establish the high speed of a type airplane 
so that it may be checked from flights of other airplanes of the 
type and under varying atmospheric conditions. The problem 
applies to an installation using a two-speed single-stage mechani- 
cal supercharger. The problem becomes more complicated when 
two-stage supercharging and exhaust turbosupercharging are 
used. Some of the factors relating to the latter type of installa- 
tion are indicated. 

The particular benefit gained from the determination of air- 
plane critical altitude is the establishment of the airplane high 
speed and the establishment of the increase of high speed effected 
by careful design of the carburetor air scoop for the attainment of 
the highest possible ram due to the dynamic pressure of the air 
entering the carburetor scoop. 


INTRODUCTION 


ITH THE ADVENT OF AIRPLANES that exceeded a 

maximum speed of 300 m.p.h., there came a 
realization that the ramming effect of air in the car- 
buretor air scoop, due to dynamic pressure, caused the 
engine b.hp. to be increased by an amount sufficient to 
increase substantially the airplane speed. 

This characteristic assumed the nature of something 
obtained for nothing and, as the value of each mile per 
hour added to the high speed is placed high in seeking 
supremacy in military aircraft design, aircraft manu- 
facturers have commenced to pay more attention to the 
effect of ram. The method to be described in this paper 
is the result of actually working out a check of the ram 
effect in increasing the altitude at which the rated horse- 
power of the engine can be obtained. 

Critical altitude will vary depending upon the power- 
operating condition of the engine—i.e., high speed in 
level flight, climb, cruising, etc. The gain is most 
significant in the high-speed operating condition. 
The present discussion is therefore limited to the critical 
altitude for high speed when maximum rated power is 
obtained from the engine in level flight. 

The problem becomes one of engine-operating char- 
acteristics to a large degree. However, since the air- 
plane manufacturer designs and installs the carburetor 
air scoop, which is responsible for the increased altitude 
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of operation, and since the actual result must be deter- 
mined by flight tests of the finished airplane, it becomes 
the responsibility of the airplane manufacturer to ob- 
tain the desired results. In the present instance, the 
need to rely on the engine manufacturer for cooperation 
in solving the combined design and flight test problem 
met with ready response from the Field Engineering 
Division of the engine manufacturer. As a result, 
this method for actually determining airplane critical 
altitude has been prepared. 


DEFINITIONS 


To clarify some of the terms characteristic of this 
problem the following definitions are given: 

(1) Ram is the increment of pressure present at the 
carburetor top deck, due to the effect of dynamic pres- 
sure of the air at the entrance to the carburetor air 
intake scoop. Ram is often measured in feet of alti- 
tude. The value so stated is the equivalent to the 
aforementioned increment of pressure under standard 
atmospheric conditions. 

(2) Critical altitude, in general, may be defined as the 
highest altitude in standard atmosphere at which the 
rated b.hp. of the engine is developed. 

(3) Engine critical altitude is the highest altitude in 
standard atmosphere at which the engine, if operated 
without carburetor temperature rise or ram, would de- 
liver its rated b.hp. at its rated r.p.m. 

(4) Airplane critical altitude is the highest altitude 
in standard atmosphere at which the engine, as in- 
stalled and as effected by ram and carburetor air tem- 
perature rise, will deliver its rated b.hp. at its rated 
r.p.m. 

(5) Supercharger compression ratio is the ratio of the 
manifold pressure to the carburetor top deck pressure. 
The manifold pressure is the pressure aft of the super- 
charger blower and before entering the cylinders. 
Carburetor top deck pressure is the pressure of the 
charge of air at the carburetor entrance. 


EFFECT OF PREVAILING ATMOSPHERIC CONDITIONS 


Engines are rated for a limited b.hp. and r.p.m. up 
to the maximum altitude at which the b.hp. can be de- 
livered. The rated b.hp. is determined by the physical 
and strength characteristics of the materials used in 
the engine consistent with satisfactory length of service. 
The highest altitude at which rated power can be ob- 
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tained is the altitude at which, because of supercharging 
by an internal or exhaust-driven blower, the necessary 
weight flow of air can be maintained to the engine. 
This altitude is called the engine critical altitude if the 
engine is operated without carburetor air temperature 
rise or ram. 

The dynamic pressure of the air entering the intake 
scoop may add to the pressure at the carburetor top 
deck. Now, since added carburetor pressure above 
normal atmospheric pressure is available at the rated 
altitude of the engine because of ram on the carburetor 
air scoop, more than rated b.hp. is available from the 
engine; or if the b.hp. is limited to its rated value, the 
full inherent pumping capacity of the engine is not 
used. This balance of the engine pumping power is 
therefore available at higher altitudes to permit ob- 
taining rated b.hp., where, without the added effect of 
ram, the rated b.hp. could not have been obtained. 

The highest altitude at which the added pressure due 
to ram is effective is airplane critical altitude, as pre- 
viously defined, and depends upon the aerodynamic 
efficiency of the air intake scoop to the carburetor. 

While ram effect of the air in the carburetor scoop 
adds to the power, carburetor air temperature rise de- 
tracts from the power, so that the net gain due to duct- 
ing the air to the carburetor is the resultant of the two 
effects. 

Airplane critical altitude differs from engine critical 
altitude because ram depends on the design of the en- 
gine cowling and carburetor air intake, and only a part 
of the total available energy of the air attributable to 
ram is available at best and a rise in carburetor air 
temperature occurs in the intake duct. If a stated 
amount of ram could definitely be obtained in any air- 
plane it would still be necessary to correct for carbu- 
retor air temperature rise. 

Airplane critical altitude is significant to an airplane 
in that its maximum speed will be obtained at the high- 
est altitude that the rated power can be obtained. 

Rated b.hp. produced at the higher altitude due to 
ram usually allows a higher speed because of the lower 
air density at the higher altitude. 

Actual airplane critical altitude can only be obtained 
from flight tests where measurements of carburetor ram 
under flight conditions will determine the efficiency 
with which the carburetor scoop has been designed for 
producing ram. 

Since standard conditions are rarely met when flight 
tests are made and since the variant factors affect the 
horsepower produced, as well as the aerodynamic drag 
and lift of the airplane, it is necessary to be able to 
correct engine b.hp. from observed atmospheric tem- 
perature and pressure conditions to standard conditions 
before the actual altitude at which rated b.hp. can be 
obtained under standard conditions, or better the criti- 
cal altitude, can be determined. Once it is known at 
what altitude rated b.hp. is obtainable, it*is possible to 
determine the high speed of the airplane with rated b.hp. 
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DESCRIPTION OF METHOD 


A method for determining airplane critical altitude 
from flight test observations follows. The determina- 
tion of critical altitude by this method is carried out by 
showing the effect of the deviations from standard 
conditions on the power calibration curve (Fig. 1). 
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ALTITUDE. - THOUSANDS ‘OF FEET 


Engine power curves corrected for ram and carburetor 
Brewster model 339, Wright R-1820-G105A. 
2, Actyal ram, uncorrected 


4, Observed data. 5, 


Fic. 1. 
air temperature. 
1, Zeroram, zero carburetor rise. 
b.hp. 3, Actual ram and carburetor rise. 
Specified b.hp. vs. altitude. 


The method as outlined here applies to a two-speed 
mechanically driven supercharged air-cooled engine on 
a single-engine tractor airplane. The observed data 
are taken from a chosen pressure altitude. Since the 
pressure altitude is observed, only temperature correc- 
tions need be made to b.hp. and intake air to correct 
to standard conditions of altitude. 

The required data are obtained by flying the air- 
plane at full throttle with best power mixture setting 
at several altitudes near the expected critical altitude, 
making b.hp. and r.p.m. measurements by means of a 
torquemeter and tachometer and observing sufficient 
atmospheric data that corrections may be made to 
standard conditions. When altitudes, corrected for the 
effect of ram and carburetor temperature rise, are 
plotted on the line of constant r.p.m. and m.a.p. (mani- 
fold absolute pressure) against b.hp. corrected for carbu- 
retor temperature rise, the intersection of the rated 
b.hp. line at rated r.p.m. with altitude on this curve is 
the engine critical altitude. In the process of develop- 
ing the corrected b.hp. curve against altitude, the ob- 
served b.hp. values are plotted against altitude. Cor- 
rections for deviations from standard conditions of tem- 
perature are made on the power curve sheet by finding 
the new b.hp. under constant conditions of r.p.m. and 
manifold pressure at altitudes represented by the 




















AIRPLANE CRITICAL ALTITUDE 


known pressure variations due to the deviations from 
standard conditions. 

Fig. 1 illustrates the differences between engine 
critical altitude, observed airplane critical altitude and 
corrected airplane critical altitude. Curve 1 shows 
the engine power corrected to standard atmosphere if 
there were no ram and no carburetor temperature rise. 
Curve 2 shows the result of correcting power curve 1 for 
the effect of ram. The correction of b.hp. for carbu- 
retor temperature rise is not yet applied on this curve, 
and the curve only represents a step in the process of 
solving the problem. Curve 3 shows the correction of 
b.hp. of curve 2 to actual carburetor temperature rise, 
and, hence, represents the true airplane critical altitude 
condition. Curve 4 shows the engine power as ob- 
tained from observed data and not corrected for devia- 
tions from standard atmosphere. Curve 5 shows the 
engine rated b.hp. against altitude as plotted from the 
engine manufacturer’s calibration of the engine. 


OUTLINE OF METHOD 


Insofar as a definite carburetor temperature rise is 
inherent in the installation, the calculation has been 
figured on the basis of this temperature rise existing 
under otherwise standard conditions. A curve of actual 
carburetor air temperature rise against outside air tem- 
perature is plotted on Fig. 2. (These data are ob- 
served data for the particular airplane installation 
being investigated.) 
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CARBURETOR TEMP RISE ~ 
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Full throttle carburetor air temperature rise. Observed 


data. 


Fic. 2. 


During a flight test for determining critical altitude, 
the following data should be observed under maximum 
rated power operating conditions in level flight: (1) 
r.p.m., (2) manifold pressure, (3) carburetor ram, (4) 
outside air temperature, (5) carburetor air temperature, 
(6) pressure altitude, and (7) torquemeter reading for 
b.hp. 

The calculations to be made from the observed data 
are as follows: (1) Determine the standard altitude pres- 
sure for the several observed sets of data by dividing 
observed manifold pressure by the supercharger com- 
pression ratio for the following conditions: (a) Ob- 
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served altitude pressure plus observed ram corrected 
to zero carburetor air temperature rise and, by dividing 
the observed manifold pressure by this corrected super- 
charger compression ratio, obtain the corrected stand- 
ard aititude pressure. (b) Corrected standard altitude 
pressure of condition 1(a) minus ram and corrected to 
include actual carburetor air temperature rise. 

(2) Determine the equivalent standard altitude rep- 
resented in the carburetor top deck for the pressures 
obtained in 1(a) and 1(b). 

Pressure of step 1(a) gives the standard altitude that 
would give the same performance as the observed per- 
formance if no ram or carburetor temperature rise 
existed for the same carburetor top deck pressure. 
These altitudes plotted against observed b.hp. would 
give the power curve for zero ram and zero carburetor 
From this curve, engine 
Pressure from the 


temperature rise condition. 
critical altitude can be obtained. 
condition of 1(b) gives standard altitudes that give the 
same performance as if ram and actual carburetor tem- 
perature rise were added. These altitudes plotted 
against b.hp. corrected for actual carburetor air tem- 
perature rise give the power curve for the airplane op- 
erating condition. 

The curve of 1(a) is plotted directly from computa- 
tions. The curve of 1(b) is obtained as follows: Cor- 
rect the altitudes of step 1(a) to standard altitudes that 
would exist when ram and actual carburetor air tem- 
perature rise are present. This is obtained by deduct- 
ing ram from the carburetor top deck pressures and 
correcting the supercharger compression ratio for ac- 
tual carburetor temperature rise over temperature of 
step 1(a) and finding the corresponding standard al- 
titudes. 

Plot these altitudes on the constant r.p.m. and mani- 
fold pressure lines of the power curves of Fig. 1 and 
draw the resulting power curve. This gives a power 
curve for actual ram and with actual carburetor tem- 
perature condition. 

The b.hp. on this last curve must now be reduced to 
correct for the effect of increase in carburetor tempera- 
ture on b.hp. itself. These corrections are made by 
correcting b.hp. to 0° carburetor temperature rise by 
the formula 


b.hp., = b.hp., V t.c.a./t.f.a. 


where 
b.hp., = corrected brake horsepower 
b.hp., = observed brake horsepower 
t.c.a. = absolute temperature of the carburetor air 
t.f.a. = absolute temperature of free air at stand- 


ard altitude conditions 


Values of carburetor temperature rise are obtained from 
Fig. 2. 

This last power curve is the curve of step 1(b). The 
altitude at which rated b.hp. occurs on this curve is the 
airplane critical altitude. This, as stated previously, 
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is the highest altitude under standard atmospheric con- 
ditions at which rated b.hp. would be obtained for the 
ram and carburetor air temperature rise conditions 
existent. 

The steps for carrying out the above procedures are 
as follows: 

(1) Note pressure altitude as read on the sensitive 
altimeter in feet. 

(2) Note the r.p.m. as set on the propeller control 
and maintained constant by governing. 

(3) Read the manifold pressure called for on the 
engine calibration curves at the selected pressure alti- 
tude. 

(4) Determine the ram, which is the measured in- 
crease in pressure in inches of water in the carburetor 
air scoop over the actual outside pressure for the oper- 
ating condition. 

(5) Note the carburetor air temperature observed on 
a thermometer in °F.—i.e., the temperature of the air 
at the carburetor top deck. 

(6) Note outside air temperature in °F. in the free 
air stream at the selected pressure altitude. 

(7) Compute b.hp. from torquemeter and tachome- 
ter readings. 

(8) Standard air temperature. 

(9) Correct b.hp. to 0° carburetor air temperature 
rise. 

(10) Convert the ram of 4 to inches of mercury. 

(11) Find the pressure that corresponds to the ob- 
served pressure altitude from the standard atmosphere 
tables. 

(12) Combine 10 and 11 to get the pressure at the 
carburetor top deck. 

(13) The compression ratio of the supercharger is 
obtained by dividing manifold pressure 3 by the car- 
buretor top deck pressure 12. The manifold pressure 
is the pressure aft of the blower and the carburetor 
inlet pressure is the pressure before the blower. 

(14) The compression ratio 13 must be corrected now 
for temperature variation of the carburetor air tem- 
perature from outside air temperature of a standard 
day. This is a function of the engine as plotted on 
Wright Aeronautical Corporation Curve No. 1693 
(Fig. 3). 

(15) Having now obtained the corrected compres- 
sion ratio as affected by deviation from carburetor air 
temperature, it is divided into the manifold pressure 3 
to determine the standard pressure equivalent to the 
conditions as they were observed. 

(16) The altitude is now found from the standard 
altitude charts corresponding to 15. This altitude is 
the altitude that under the observed conditions would 
have been obtained if there were no ram and no de- 
viations from standard pressure. 

(17) Plot the altitudes 16 on constant r.p.m.-m.a.p. 
line from b.hp. 9. The intersection of this curve with 
the rated b.hp. of the engine gives the engine critical 
altitude. 





U | 
20-10 O 10 20 30 40 50 60 70 80 9 WoO 
CARBURETOR AIR TEMPERATURE ~ °F 


Fic. 3. The ratio of pressure at blower outlet to pressure at 
impeller inlet as a function of carburetor air temperature is P2/P;. 
W.A.C. curve No. 1693. Formulas: P2/P,; = (y + 1)%-53, y= 
v?m/5960TR, where n = 0.67. 


(18) Subtract the ram in inches of Hg from the 
standard altitude pressure 15. , 

(19) Determine the standard altitude represented by 
the standard pressures 18. 

(20) Find the carburetor air temperature °F. at the 
altitude 19 from Fig. 2. 

(21) Correct the compression ratio 14 for the carbu- 
retor air temperature 20. 

(22) Find the standard altitude pressure for the cor- 
rected compression ratio 21 by dividing the pressure 
ratio into the manifold pressure 3. 

(23) Subtract the pressure due to ram from the pres- 
sure 22. This gives the pressure representing the stand- 
ard altitude that would prevail if there were ram and 
actual carburetor temperature rise. 

(24) Find the altitudes obtained from pressures 23 
and plot on lines for variation of b.hp. versus altitude 
as constant r.p.m. and manifold pressure on Fig. 1. 
This curve represents the b.hp. uncorrected for car- 
buretor temperature rise at the correct altitudes. 

25) List b.hp. at altitudes 24 from Curve 2 on 
Fig. 1. 

(26) List standard temperature of altitudes 24. 

(27) List carburetor air temperatures from Fig. 2 for 
altitudes 24 with temperatures in absolute degrees. 

(28) Correct b.hp., = b.hp., V t.c.a./t.f.a. 

(29) List corrected b.hp. and plot against altitudes 
24. 

(30) Where rated b.hp. intersects curve 29 is the air- 
plane critical altitude. 
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TABLE 1 


Sample Calculations 








Test 2 Test 3 
1. Observed pressure altitude—ft. 21,350 19,700 18,800 18,500 
2. Engine—tr.p.m. 2,355 2,355 2,355 2,355 
3. Absolute manifold pressure, P, in. Hg 30.7 33.0 35.0 35.3 
4. Ram—in. H,O 14.0 17.0 19.3 19.3 
5. Carburetor air temperature, 7;—°F. 19.5 23.0 25.0 30.0 
6. Free air temperature, 7,—°F. —13.0 —4.0 +1.4 +3.0 
7. Observed b.hp. 696 753 808 810 
8. Standard air temperature, 7,—°F. -17.1 11.2 8.0 —7.0 
9. C.hp. (b.hp. corrected from 7; to 7, 724.5 781 837 842 
10. Ram—in. Hg 1.03 25 1.42 1.42 
11. Pressure altitude—in. Hg 12.98 13.92 14.45 14.63 
12. Carburetor top deck pressure, P;—in. Hg 14.01 15.17 15.87 16.05 
13. Observed pressure rise ratio—(P2/P;) obs. 2.19 2.18 2.21 2.20 
14. (P:/P) std. (corrected from T; to 7,) 2.265 2.25 2.28 2.28 
15. =. = (standard alt. press. for zero ram and 7, — 7, rise) 13.55 14.67 15.35 15.48 
(P2/P;) std. 
16. Altitude (from 15)—ft. 20,350 18,430 17,330 17,120 
17. Plot c.hp. from 9 against alt. of 1 and follow constant m.a.p.-r.p.m. line to altitude 
from 16. See Fig. 1. Final c.hp. is: 716 771 825 830 
From Fig. 1 the engine critical altitude is 17,800 ft. 
18. Standard altitude pressure from 15 minus ram from 10—in. Hg 12.52 13.42 13.93 14.06 
19. Pressure altitude corresponding to 18 5 22,200 20,570 197 194 
Standard air temp., 7,, for 19—°F. ; ~20.0 14.4 —11 ~10 
20. Carburetor air temperature, 7) for 19 from ‘Fig. 2—°F. 17.0 19.6 a 21 
21. P2/P;[(P:/P;) std. from 15 corrected to T; of 20 from T, of 19] 2.195 2.185 2.22 °: 
> 
22. Carburetor top deck pressure, Pi, = =. ae 14.01 15.10 15.76 15.90 
(P2/P;) of 21 
23. FP; from 22 minus ram from 10 = altitude pressure 12.98 13.85 14.34 14.48 
24. Altitude (from 23)—ft. 21,360 19,816 18,984 18,750 
25. Hp. at altitude of 24 (from final c.hp. of 17 along constant m.a.p.-r.p.m. line to 
altitude of 26), Fig. 1 725 782 838 844 
26. Standard air temperature, 7, for 26—°F. 7.2 -11.66 —8.64 —7.9 
27. Carburetor air temperature, 7, for 26 from Fig. 2—°F. +18 +21 +19.4 +19.1 
28. Hp. of 25 corrected from 7, of 26 to 7; of 27 698 762 815 820 
29. Plot hp. of 28 against altitude of 24 on Fig. 1 
30. From Fig. 1 the airplane critical altitude is 19,200 ft. 
SAMPLE CALCULATION CONCLUSION 


A sample calculation is presented from data observed 
on a single-place single-engine fighter with nine-cylinder 
air-cooled radial engine (Table 1). The air is taken 
into the carburetor by a ramming scoop at the leading 
edge of the nose ring. The entrance to the scoop is 
faired to permit a minimum of turbulence at the inlet. 
The ram represents an increase of 17 in. HO pressure 
at the carburetor top deck as against a dynamic pres- 
sure of 25.3 in. before the entrance. Hence, the scoop 
effectiveness is 67 per cent of the total available pres- 
sure increase neglecting slipstream. 

The increase in altitude for rated b.hp. is 2,100 ft. 
above that rated for the engine and 1,400 ft. above the 
actual engine critical altitude. 

The corrected high speed of the airplane is 327.5 
m.p.h. This is 7 m.p.h. in excess of what it would be at 
rated altitude of the engine, 3 m.p.h. of which is due to 
the engine exceeding its rated altitude and 4 m.p.h. due 
to ram on the carburetor intake. 





This method of determining airplane critical altitude 
is not difficult to apply when the curves of effect of 
carburetor air temperature versus compression ratio 
are available. The method as used did not require 
special flight test equipment installed other than nor- 
mal equipment for routine flight tests. It has been 
found possible to check several airplanes of this same 
type on different occasions to give the same ram effect 
and consistent high speeds. A later model of the air- 
plane was quite satisfactorily checked against esti- 
mated ram from wind-tunnel tests. 


FUTURE DEVELOPMENT 


It is unfortunate that this problem as presented 
applies to an installation that is not in line with the 
latest developments where its application is important. 
The development of two-stage and exhaust turbine- 
driven superchargers on high-speed aircraft presents a 








Test 4 
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modification of the problem, the development of which 
has not yet been simplified. 


In the above-mentioned installations an auxiliary 
stage of supercharging is involved, the corrections for 
which have to be made before the carburetor top deck 
is reached by the air charge. Beyond the carburetor 
top deck in the system, the problem remains essentially 
the same. 


There is a certain pressure drop that takes place 
through the auxiliary stage and intercoolers in the two- 
stage and turbo systems which depends on the rate of 
flow through the system and the efficiency of the 
intercoolers involved. 


When curves are available from which the compres- 
sion ratio of the auxiliary supercharger are established 
and the pressure drop across the intercooler installation 
can be established, a method can be worked out 


whereby, with the aid of additional observed pressure 
and temperature data, corrections can be obtained 
which will permit determining the actual critical alti- 
tude of these installations in an airplane. 

It will be desirable to standardize the temperature 
at the carburetor top deck to which the air through the 
intercoolers must be cooled in order to simplify the 
method. At present, different temperatures are al- 
lowed at the carburetor top deck by different agencies 
working on the development of the system, and as- 
sumed values of pressure loss are used for the system 
between the air intake and the carburetor top deck. 

At best a solution of this problem will be more com- 
plicated, probably commensurate with the added com- 
plications of the supercharger installation itself. How- 
ever, when more data have been established on the 
supercharging systems, a workable method of deter- 
mining airplane critical altitude will be developed. 





Erratum 


In the “Letter to the Editor’ by G. A. Brewer 
(JOURNAL OF THE AERONAUTICAL SCIENCES, August, 
1942, pp. 371 and 372), Fig. 1 was lost in transit and 
was omitted. This figure is published herewith. 
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LELMENT NQ / LLEMENT MUTUIMLLY = Of ACA 


Fie. 1. 


On page 371 of the August JouRNAL the illustration 
shown as Fig. 1 should be Fig. 2a, and on page 372 


Fig. 2 should be Fig. 2b. 


In the second paragraph, 


second column, page 372, the second sentence should 
be changed to read: “‘Referring to Figs. 2a and 2b, 
the nonuniform nature of strain distribution along a 
specimen of polycrystalline material is clearly shown 
in the outlined zone.” 
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Vector Solution of the Three-Degree Case 
of Wing Bending, Wing Torsion, Aileron 


Flutter 


LEE ARNOLD* 


Republic Aviation Corporation 


INTRODUCTION 


ie HAS BEEN SEEN by investigators that the two- 
degree analysis based on Theodorsen’s!: * theory is 
inadequate for the accurate determination of the 
critical flutter speeds of an oscillating airfoil-flap sys- 
tem. Various graphic methods of solution have been 
devised to eliminate the computations required in the 
solutions of the various two-degree cases,*: * but until 
now no simplified solution has been developed for the 
Though the method of solution of 
the three-degree case is straightforward and routine in 
character, the actual inherent complexity of the com- 
putations involved is such that it warrants the following 
simplified vector method of solution. 

The solution is based on the fact that the components 
of the vibratory determinant are complex in form (real 
part plus an imaginary part) and therefore can be repre- 
sented by a vector in the complex plane. For a given 
set of physical parameters the terms of the determinant 
are functions of the auxiliary flutter parameter k = 
wb/V and w (w is the flutter frequency, } the airfoil 
semichord, and V the flutter speed), while the remain- 
ing terms are functions only of k. Expansion of the 
determinant results in terms that are functions of k 
and w and terms that are functions only of k. Since 
both of these groups are complex in form, they can 
still be represented by vectors. 

Transposing these vectors, one on either side of the 
equality sign, the flutter speeds will be found by deter- 
mining those values of k and w for which the two vectors 
are equal. 


three-degree case. 


SYMBOLS 
b = half chord of airfoil at reference station 
be = distance of flap hinge aft of mid-chord 


ba = axis of rotation aft of mid-chord 


p = mass of air per unit volume 
M = mass of wing per unit length 
Sa, Sg = static moments of airfoil per unit length of airfoil 


flap about a and flap about c 
Ia, Ig = moments of inertia per unit length of airfoil flap 
and flap about a@ and ¢, respectively 
apb?/M = the ratio of the mass of a cylinder of air 
of diameter equal to the chord of the airfoil to the 


Presented at the Aerodynamics Session, Tenth Annual Meet- 
ing, I.Ae.S., New York, January 30, 1942. 

* Now associated with Pratt & Whitney Aircraft Division, 
United Aircraft Corporation. 
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mass of the airfoil, both taken for equal length 
along the span 


bra = V/ Te M = the radius of gyration 

bxe = the center of gravity distance of the airfoil from a 
Wa = the frequency of torsional vibration about a 

brg = +/Ie/M = reduced radius of gyration of the flap 
bx = Sg/M = reduced center of gravity of flap from « 
wp = frequency of torsional vibration of flap about c 
Wp, = frequency of airfoil deflection 

T; = functions of c 

V = critical flutter speed 

w = critical flutter frequency 

k = wh/V 


MATHEMATICAL DEVELOPMENT 


The equations of motion resulting from the combi- 
nation of the inertia, elastic, and aerodynamic forces 
and moments of an oscillating airfoil-flap system, Fig. 
1, at the state of steady oscillation which separates 
the condition of stability and instability are: 


(Raw + ilaa)ao + (Rap + ilap) Boe + ) 
(Ran + tIqn)hoe' = 0 

(Roa + tT a) a0 + (Rog + iL yp) Boe + 
(Ron + tIon)hoe™ = 0 

(Rew + 4 co.) 0 -+- (Reg + iI.g) Boe + 
(Ren + iT »)hoe™™ = 0 


| 
La) 


where a, 8, and fp are the maximum amplitudes of the 
airfoil about the torsional axis, the flap about its hinge, 
and the vertical displacement of the system, respec- 
tively. ind 
Since the above equations are homogeneous, in order 
that solutions other than the trivial zero solution exist, 
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it is necessary that the determinant of the coefficients 
of Eq. (1) vanishes, so that 
| - > tl oa Rag a 1 ag Ran + tT qn| 
Rra + tDra Rog + iT yp Run + Tn) = 9 (2) 
Rea 43 1 Reg i il eg Rx + il en 


The components of the determinant above are of the 


form 


‘dD *F Ya" . To* Wa” , 
(Row + tl en) oe = oan + 1 La) ~~ — Ros + 
K w* 
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where the (Rny’)’s and (Im,)’s are of the form 


} my = (1/3 + a’) + 2(1/4 sii a*)G/k ia 
2(7/2 + a) F/k? 
- 
Rap’ = — (T1+ (€ — a)T, + (14+ Ti) (1/k*) + 
ri 
(a + '/2)[(TuG/k) — (2TwF/k?)]} 
Ray’ = a + 2(a + '/2)G/k 
1 
Rea’ = — {T7 + (¢ — a)T, + Trl F/R? — 
vs 


(*/2 — a)G/k]} 


1 
Rye’ = AE + (Ts — 147 )1/k? — 
a 


T 
ne — G/k — T»F/k? )] 


Ron’ = 1, w[T) ” T2G/k] 





Re’ =a — “dea — a)G/k + 2F f (4) 
Rep’ = 1/x(T, — (TunG/k) + 2T Fe 
Ra’ = — [1+ (0G ‘)) 
Toe = — 1/R{2(a + 1/2) [2/2 — a)F + G/k] — 

("/2 — a)} 
Iga = — 1/mki(a + 3/2)(TuF + 2TwG/k) + 

2p + (1/2 — a)T,} | 

Ign = — 2(a+ '/2)F/k 
Tha = 1/k\ Ty[('/2 — a) F + G/k] +p —Ti — 


7/2) 
Ing. = 1/0? [(TuT 2/2) F/k + TiT2G/k?] 
Ty, = T12/4F/k 
Tea = 1/k[2(*/2 — a) F + 2G/k + 1] 
Teg = 1/mk[TuF + (2TwG/k) — T;) 
in = ZF 








kK 
a 
(Rap + thas) = - tle ~ + Rap’ +ilap | 
K K | 
ie | 
(Ran + tT qn) = a + R,,’ + 7 | 
kK } 
NS ee: ae i 
K K 
, a | 
" ms rg” rg” We" 
Bue + the = — — + 11 +800 a er + Roe’ + (3) 
. 2 
tIyg | 
, xp , : 
Ron + tly, = — — + Ry’ + tly, 
K 
Xa . 
Re + rt = es —-+ a + tL cx 
K 
7 vp P ; 
Rep + thep = — a Reg’ + t1ee | 
| 
a : 
Ra + ila = — + ee ign) - = E+ Be! + iy, J 
Making the substitutions 
01 = -9 a (c — a) = | 
Kk | 
ee 2 (1 -% _ . (5) 
K o) K | 
l ( =) Xa 
os = -\1 — laa } 
K Ww K 
into Eq. (3), there results 
e ™ Vos 
Res + tl on = igal or i ) » 
K 
(Rea’ + iTea) | 
- - ' rp’ \ | 
Rop + tle = — o2 — 1gg\ o2 — 3 = 9 
(Rop’ + ilve) | 
= = 1 
Rear la = — —1 ( — *) er 
at tle G3 12h a + 4 (6) 
(Ren’ + tI en) 
Rap + tlap = — m + (Rap’ + tag) 
Ria + tla = — 1 + (Rea’ + tlre) | 
Ro + thy» = — m2 + (Ron’ + iI yn) 
Reg + il ep = oa + (Rep’ + 11.) 
Rar + tI an a + (Rar’ + tI qn) 





Rea se OT co a ae + ea + 1D cx) ) 


For fixed values of c and a, the above R,,,,’ and Im 
terms are functions only of k, so that each term Ry» + 
tImn is made up of a term containing the inertia, elastic 
and friction effects, and a term Ry’ + iIm_, dependent 
only upon the auxiliary flutter parameter, k. 

Since each term of Eq. (6) is complex in form, it 
can be represented by a vector in the complex Ryn, 
Imn plane. These vectors will be of magnitude A,,, and 
will form the angle Wm, with the R,,., axis. 

Each of these vectors will take the form 


‘ iym? 
) an + 1 _ has” 


and Rin + tlen = Anati™ (7) 


Substituting Eq. (7) in Eq. (2), expanding, and then 
placing all barred terms on the right-hand side and un- 
barred on the left, there results 


AgpAmA cae Veet eetven) 4 Ay A pA ane ra treat ¥en) 
AcAedgdOuctontin 4 Arededgintvatia) 4 
A A apArae’ Yt a8 tH) 4. A AypA ene’ Vert pt ven) 

i aa f1¢ 


This has the form 
Hie + Hee™ = Aye™ ae Aye oo Ae™ +- Free’ 
or He® = He* (8) 
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WING BENDING, WING 


CONSTRUCTION AND USE OF VECTOR CURVES 


Since the procedures for the construction of each of 
the vector plots Rnn + iIm, are similar, only the plot of 
Roa + ica Will be explained completely. It is recalled 
that 
7 - Ta” ' 
Roa + thon i oa ital o1 as ) + (Rea’ + tI oa) 

\ K 

Since the expression Rye’ + ila is purely a function 
of k for a fixed c and a, it can be plotted independently 
of the physical parameters as shown in Fig. 2. 
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The curves Iga vs. Rya’ and Iga vs. 1/k are plotted 
for permanent use. (These curves have been plotted 
fora = —0.1, —0.2, —0.3, —0.4, —0.5.) A straight 
line is drawn through the points, A measured a positive 
distance r,”/x on the real axis, and B measured a dis- 
tance —gar.’/k on the imaginary axis. A point is 
dropped from Rye’ = o; in a direction perpendicular to 
the R,.’ axis to the line AB to form the origin, 0, of the 
vector. Then from the chosen value of 1/k on the 1/k 
axis a point is determined on the Iga vs. Raa’ curve 
which when transmitted parallel to the R,.’ axis will 
determine the point P, the end point of the vector, on 
the Iga vs. Rea’ curve. 

Identical procedure was followed for the determina- 
tion of the vector 


Rye 4. ilys s=-—- & = iga(oe —- r”/k) + (Ryp’ + ilyg) 


The point A will be at (rg?/x, 0), the point B at 
(0, — gere*/x), and o2 will be used instead of o. 
Also for the determination of the vector 


Re + il en = 1gn(os _ (1/x)] + (Ren’ + iT x) 


The point A will be at (1/x, 0), the point B at (0, 
—g,(1/x)), and o3 will be used instead of o;. 01, o 


TORSION, 
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and o3 have a definite relation to each other and cannot 
be chosen at random. 

For a given value of o, the corresponding oz and a; 
are determined from the system shown in Fig. 3. 


a a 
T Y7 
A (#2) 
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For the unbarred terms 
Ros + | — “a + § ay + tI mn) 


Run’ + tImn has been plotted for a complete range of 
k’s for various c and a combinations. The procedure is 
similar to that employed in determining the barred 
vectors except that the vector origin is located at the 


point Ring’ = ni 


PROCEDURE 


1. The lengths A,,, and the angles y¥,,, are measured 
and tabulated for a few values of k in a practical range. 

2. The lengths Amn and the angles Ym», are measured 
for a range of o;’s (and the corresponding o2’s and o;3’s) 
and a range of k’s to be used with each value of o;. 

3. The values of Amn, Wmn are substituted into the 
unbarred terms of the left side of Eq. (8) forming the 
vector He”, and then the values Amn, Ymny Amn 200 Yn 
are substituted into the right-hand side of Eq. (8) 
forming the various vectors He”. 

4. Plot H vs. @ for each value of k, and H vs. k as 
shown in Fig. 4. 

5. Plot Avs. 6 fora fixed k = k and varying o;’s, and 
6 vs. a1. 

6. Superimpose the H, @ curve onto the 7, 6 curve. 

7. Translate P;, P2, P3, ..., P, horizontally to the 
H vs. k curve and then down to the ordinate k;, Re, ..., 
k,, respectively, in the k, k plane. These (k, k) points 
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will be the loci of a curve in the k, k plane. Where the 
bisector of the & and & axes intersects the k, k 
curve, the point M is determined. 


8. Translate P;, Po, ..., P, vertically downward 
to the k curves of the o; vs. 6 group. These points on 
ky, ko, ..., Ry will form the locus of a (k, 1) curve. 


9. Translate M vertically upward to the point N 
on the H/ vs. k curve, then from N horizontally to R 
on the H vs. @ curve, then down vertically to S on the 
k vs. o; curve, and finally horizontally to T on the o; 
axis. 


—> 
10. Where MN intersects the & axis k is determined, 
and the point T determines o. 


11. From the relations a, = (fq?/x){1 — (we?/w?)| 
and k = wb/V, there results 
5 Wad l 
Y= —- = (9 


k Vi — [oy CA 9) 


from which the flutter speeds are determined. 
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Predicting Stiffener and Stiffened Panel 
Crippling Stresses 


HAROLD B. CROCKETT* 
Lockheed Aircraft Corporation 


SUMMARY 


An empirical method is presented for predicting the crippling 
stress at L/p < 20 of 24S-T Alclad bent-up sheet stiffeners, 
24S-RT Alclad bent-up sheet stiffeners, and 24S-T extruded 
stiffeners acting alone or in conjunction with 24S-T Alclad sheet 
asa panel. It is also shown that predicted column curves may 
be constructed by drawing a straight line tangent to the Euler 
column curve (c = 1.5) through the predicted crippling stress at 
L/p = 20 (except through L/p = 40 for extrusion stiffened 
panels); the predicted stress at L/p = 20 (or 40) locates the 
horizontal cutoff for the section of the column curve with L/p < 
20 (or 40). 

The general method, which consists of summing the crippling 
loads of the component elements of the stiffener and then divid- 
ing by the summation of the elemental areas to obtain the crip- 
pling stress, is discussed in several texts on airplane structural 
analysis; however, enough test data have not been available 
heretofore to establish general empirical curves for predicting the 
strength of any stiffener or stiffened panel of reasonable propor- 
tions. 


INTRODUCTION 


a STRUCTURAL ENGINEERS continually are 
faced with the problem of finding the crippling 
(local compressive failure) strength of aluminum alloy 
stiffeners and stiffened panels. If test results are not 
available, the engineer ordinarily estimates the total 
load the stiffener will carry by applying theoretic equa- 
tions (such as those in references 1 and 2) to determine 
the buckling loads of the various flat and curved ele- 
ments of the stiffener and then summing the loads car- 
ried by the elements. This paper gives an empirical 
method for predicting the crippling stress of stiffeners 
and stiffened panels by using design curves that are 
trial-and-error variations from the theoretic to give the 
best possible agreement with test data. 

It is suggested that one read this paper entirely 
through before attempting to use the curves and tables 
given in the first part. Especially should all limitations 
of L/p, b/t, etc., be studied, as well as the notes ap- 
pended to Table 1. Worth-while emphasis may be 
placed on the fact that Table 1 is a compact device for 
showing cross-section sketches /o scale of all the speci- 
mens offered as the basis for this report, besides demon- 
strating typical examples of stiffeners to which are ap- 
plicable the “stability shape factors.’’ These factors 
are general empirical coefficients indicating the rela- 
tive strength of amy (1) vertical-legged and (2) slanting- 
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LS 219, 208 LS 209, 200,210 15236, 232 L5264 
NOoTEs: 
*(a) The stiffeners shown are those tested. K applies to any 


stiffeners judged to lie within the range of the method given. 
(b) One division represents '/,in. (c) The underlined identifi- 
cation numbers indicate that those stiffeners were tested with, 
and in some cases also without, sheet attached. 

* K is given for stiffeners of various materials and angu- 
i of legs assuming that thickness of sheet is less than, or 
equal to, thickness of stiffener unless otherwise indicated. 


legged stiffeners judged by the reader to lie within the 
range of the curves given. 

Many engineers will not care to extrapolate the design 
curves of Figs. 2 and 3 nearly as far as shown; how- 
ever, in this respect, the following preliminary state- 
ments may be of assistance: 

1. The extrapolated portions of the design curves 
for the flat elements of stiffeners (Fig. 2) fair into the 
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(b) 


Fic. 1. Conventions. (a) Sheet and extruded stiffener con- 


panels assumed simply supported (no edge free) are 
increasingly conservative for b/t > 45 and that at b/t 
of 100 these stiffened panels act as would theoretically 
fixed-edge panels. Flat elements of stiffeners act 
similarly. 

2. The extrapolated portions of the design curves 
for the no-edge-free curved elements of stiffeners (Fig. 3) 
lie completely within the range of, and indicate compres- 
sive strength values only a little greater than, those 
curves for circular corrugations given in Fig. 5-4 of 
reference 1. In fact, good agreement exists between the 
column stress at L/p = 20 given by Fig. 5-4 of reference 
1 for circular corrugations and the crippling stress 
computed using the present design curves and the logical 
assumption that one-half pitch-width of a circular cor- 
rugation acts as two circular quadrants on either side 
of a flat element two quadrants wide (all three elements 





ventions. The convention for determining the areas of elements : . : . 
. ° ° ° : ° ° Ss ) ave yo s1qaes § TS O eGo y 
of extruded stiffeners is to be used in connection with predicting assumed to have two sides simply “=r rted, i.e., with 


the average crippling stress; the load carried by the stiffener is no edge free). 


the product of the average stress and the total area of the 


stiffener. (b) Sketch for special rule given. 
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Fic. 2. Crippling stress of flat elements of 24S stiffeners in 


compression. 
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Fic. 3. Crippling stress of circular elements of 24S Alclad 


sheet stiffeners in compression. 


3. The second part of the paper contains the com- 
plete test results upon which the curves of Figs. 2 and 
3 are based. The quantity of data presented should 
instill the engineer with confidence in the design curves 
and afford him the opportunity to correlate his own 
private method for predicting stiffener strengths with 
these tests. 


PREDICTING CRIPPLING STRESSES: L/p < 20 


General Method 

The crippling stress of a stiffener may be predicted by 
summing the crippling loads of the component elements 
of the stiffener and then dividing by the summation of 
the elemental areas, as expressed by Eq. (1): 


D Data Sen 


> Dale 


F.. = K = KF. (1) 


where 


F., = final predicted crippling stress, lbs. per sq.in., 
for L/p < 20 


F..; = uncorrected predicted stress = 2 Data fen/Z Ontn 
K = the stability shape factor given in Table 1 
batn = area of individual element, sq.in. 

fon = average ultimate stress of the individual ele- 
ment, given empirically by Figs. 2 and 3 for 
flat and curved elements, respectively, Ibs. 
per sq.in. 

p = the radius of gyration of the stiffener alone 
about an axis parallel to the sheet in a stiff- 
ener-sheet combination, in. 

L = length of stiffener or panel, in. 


The crippling stress of a stiffened panel may be deter- 
mined as above by considering the effective width of 
sheet acting with the stiffener to be a simply supported 





classic theoretic buckling curves as }/t approaches 80. 
Results of Lockheed tests detailed by Fischel* show that 
theoretic curves giving the strength in compression of 


that is, no-edge-free—element of the stiffener. The 
reader will recall that by definition the effective width 
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Fic. 4. Effective widths at various stiffener stresses of 24S 
Alclad sheet incompression. Key: (1) 24S-RT Alclad sheet; (2) 
24S-T Alclad sheet; (3) 24S-RT Alclad sheet and stiffeners; (4) 
24S-T Alclad sheet, 24S-RT Alclad stiffeners; (5) 24S-T Alclad 
sheet and stiffeners; (6) 24S-T Alclad sheet, 24S-T extruded 
stiffeners. 


of sheet is said to act at the average stiffener stress. 
This is the basis upon which the effective width is esti- 
mated from Fig. 4. However, to predict the panel 
crippling stress the effective width is then assumed to be a 
no-edge-free element of the stiffener-sheet combination and 
to act at the much lower than average stress thereof. 
The procedure is discussed fully in another sec- 
tion. 

The conventions adopted for measuring 3, b,, t, and 
R are shown in Fig. la, where 


b = developed width of element, in. 

b,, = stiffener spacing (width sheet), in. 

b, = effective width of sheet acting with stiffener, 
in.; estimated from Fig. 4 

ts; = thickness of stiffener, in. 

t,, = thickness of sheet, in. 

R = the centerline radius of a curved element, in. 


Fig. la also indicates clearly the significance of the 
terms ‘‘one-edge-free’’ (one edge simply supported) 
and ‘‘no-edge-free”’ (two edges simply supported) stiff- 
ener elements. 

An arbitrary rule is stated as follows to take care of 
instances when a narrow flat element lies between two 
curved elements: If two curved elements of equal 
radius (centers on same side of sheet) are separated as 
shown in Fig. lb by a flat element with width b < R, 
the average stress of the three elements is taken as that 
of a curved element of radius R’ = R + (6/2). 











CRIPPLING 


STRESSES 


Predicting Crippling Stress of Stiffener Acting Alone 


It is desirable first to give a detailed example of pre- 
dicting the crippling stress of a stiffener acting alone 
and then to consider a stiffener with sheet attached, be- 
cause the techniques are slightly different. 

An outline of the steps involved in predicting the 
strength of a stiffener acting alone follows: 

(1) Make a table as shown in Fig. 5. (a) Enter 
sketch, 5, nominal R, actual R/t, b/t. (b) Obtain f,, 
from Figs. 2 and 3. (c) Enter Data, Ontn fon, Z Ontn, and 


D Dntnfen. (A) Compute Fyog = 2 Data fen/= Ont. 
(2) Find the stability shape factor in Table 1 and 


compute F,, = K Fic}. 
Example 1: Given 0.040 gauge LS 160 24S-T Alclad 





section. Find the crippling stress of a column with 
L/p = 20. 
(1) Makea table as shown in Fig. 5. 
z Data Sen 2,311 ’ 
Fee = = zal = 50807” 28,600 Ibs. per sq.in. 


From Table 1, for a vertical-legged 24S-T Alclad stiff- 
ener, K = 1.07. 
(2) Therefore, the crippling stress 


Fee = KFeex = 1.07 X 28,600 = 30,600 Ibs. per sq.in. 


For this stiffener the corrected test crippling stress 
F.,' = 32,400 Ibs per sq.in. (see Table 2), or the test 
stress is 6.0 per cent of the predicted stress higher than 
the predicted stress. 


Predicting Crippling Stress of Stiffener Acting with Sheet 


Ordinarily the addition of sheet to a stiffener in com- 
pression reduces the crippling stress of the stiffener be- 
cause, as the sheet buckles, wrinkles tend to form across 
the stiffeners, even though prevented from actually 
doing so by proper rivet or spot-weld spacing. This 
tendency causes a side load on the stiffener which re- 
duces the crippling stress. By assuming that the 
effective width of sheet acting with the stiffener becomes 
a no-edge-free element in the stiffener-sheet combina- 
tion, it is possible to estimate the crippling stress of the 
combination. The effective width to use as a no-edge- 
free element may be determined approximately from 
Fig. 4, the development of which is discussed by Fischel,* 
by assuming the stress of the combination to be F,,; of 
the stiffener alone. The latter assumption is accurate 


7 8 9 


1 2 3 4 5 6 11 
Part R’ R/t 
No. n - b, ba/ta Nominal Actual Sen bds Datnfen Fret 
1 0.040 0.646 3/16 5.2 28,200 0.0258 729 
; 2 0.040 0.648 16.2 33,700 0.0259 74 
3 0.040 0.225 1/s 3.6 43,600 0.0090 392 
4 0.040 0.500 13.5 15,800 0.0200 316 











28,600 
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TABLE 2 


24S-T Alclad Stiffeners, 24S-T Alclad Sheet 











1 2 3 4 5 6 7 8 
Fc’, Lbs. Predicted Margin, 
Stiffener Sheet f per Sq.in., Test F..1, Lbs. per Fee = K Fei, Cx ~ 1)100 
Stiffener* Gauge, ts: Gauge, ter Corrected Sq.in. K Lbs. per Sq.in. 3 ) ; 
LS 160 0.025 23,700 21,400 1.07 22,500 5.3 
0.032 0.020 28,300 24,100 1.07 25,800 9.7 
0.032 0.025 27,500 23,400 1.07 25,000 10.0 
0.032 0.032 25,800 22,400 1.07 24,000 7.5 
0.032 0.040 23,500 21,300 1.07 22,800 3.1 
0.040 32,400 28,600 1.07 30,600 6.0 
0.040 0.020 30,300 26,400 1.07 28,200 7.4 
0.040 0.025 28,900 26,100 1.07 27,900 3.7 
0.040 0.032 27,500 24,400 1.07 26,100 5.4 
0.051 0.020 36,500 30,200 1.07 32,300 13.0 
0.051 0.025 34,500 29,600 1.07 31,700 8.8 
0.051 0.032 32,200 28,400 1.07 30,400 6.0 
0.051 0.040 33,100 27,100 1.07 29,000 14.1 
0.064 43,300 35,000 1.07 37,400 15.8 
0.064 0.025 37,300 32,800 1.07 35,100 6.3 
0.064 0.032 32,300 31,600 1.07 33,800 — 4.4 
0.064 0.040 32,600 30,400 1.07 32,500 0.3 
0.081 45,300 38,400 1.07 41,100 10.2 , 
W = */,as above 0.032 0.040 24,050 21,300 1.07 22,800 5.5 
W = 3/, 0.032 0.040 22,000 19,700 1.07 21,100 4.3 
We=1 0.032 0.040 17,400 19,100 1.07 20,400 —14.7f 
LS 161 0.025 21,500 20,400 1.07 21,800 — 1.4 
0.040 29,300 26,900 1.07 28,800 .8 
0.064 41,600 34,500 1.07 36,900 12.7 | 
0.081 45,500 38,200 1.07 40,900 11.2 
LS 166 0.051 0.032 35,500 32,100 1.07 34,400 3.2 
XLS 165 0.025 20,700 19,100 1.07 20,400 1.5 
0.040 26,600 25,100 1.07 26,800 — 0.8 
X-785-A 0.025 19,600 17,200 1.07 18,300 Jen 
0.040 27,000 25,400 1.07 27,200 — 0.7 
0.064 41,100 36,800 1.07 39,400 4.3 ; 
0.081 46,000 42,000 1.07 45,000 2.2 
X-785-B 0.025 18,400 15,600 1.07 16,700 10.2 
0.040 25,600 23,200 1.07 24,800 3.2 
0.064 34,900 32,600 1.07 34,900 0 
0.081 44,200 38,700 1.07 41,400 6.8 
X-785-C 0.025 19,600 17,500 1.00 17,500 12.0 
0.040 26,200 26,200 1.00 26,200 0 
0.064 37,700 37,400 1.00 37,400 0.8 
0.081 45,200 41,200 1.00 41,200 9.7 
X-785-D 0.025 22,000 20,500 1.07 21,900 0.5 
0.040 31,100 29,100 1.07 31,100 0 
0.064 42,900 39,600 1.07 42,400 1.2 ' 
0.081 : 47,500 42,400 1.07 45,400 4.6 
LS 140 0.040 29,000 27,400 1.07 29,300 — 1.0 
LS 173 0.032 0.025 26,500 25,600 1.00 25,600 3.5 
0.032 0.032 26,100 24,300 1.00 24,300 7.4 
LS 174 0.032 0.025 27,400 26,400 1.00 26,400 3.8 y 
0.040 0.040 19,600 25,200 1.00 25,200 —22.2** t 
0.051 0.032 32,400 28,500 1.00 28,500 13.7 s 
0.051 0.040 32,400 27,900 1.00 27,900 16.1 
0.064 0.032 37,100 31,900 1.00 31,900 16.3 s 
* Sketches shown drawn to scale in Table 1. 
+t No sheet gauge given indicates stiffener tested alone. S 
t Stiffener of poor proportions. , 


** Limiting thickness in Table 1 based on this point. Test might be bad, although three points on column curve agree well. 











PREDICTING CRIPPLING 


enough considering the empirical nature of the method. 
Since 0,/t;, ordinarily is about 40, the stress of the ef- 
fective element is low. Asa result, the predicted aver- 
age crippling stress of the stiffener-sheet combination is 
lower than that of the stiffener acting alone. 

An outline of the steps involved in predicting the 
strength of a stiffener acting with sheet follows: 

(1) Estimate F,,; of the stiffener alone as detailed 
above. 

(2) Knowing the dimensions of the panel, use Fig. 4 
to find the effective width of sheet b, acting with the 
stiffener at the stress F,,;. 

(3) Calculate b,/t;, and determine f,, for the no- 
edge-free effective sheet element. 

(4) Calculate bjt, and Ontrfon for the effective sheet 
element; add these to the = d,t, and 2 Opty fn, respectively, 
obtained in step 1 (table similar to that of Fig. 5). 

(5) Compute F,,; of the combination 


New Feet =2 Datnfen/Z Datn 


(6) Find the stability shape factor in Table 1 and 
compute F,, = KF ci. 

Example 2: Given the 0.040 gauge LS 160 24S-T 
Alclad section of Example 1 acting with 0.025 gauge 
24S-T Alclad sheet in a panel having ),,/t,, = 94 and 
L/p = 20. Find the crippling stress of the stiffener, 
effective-sheet combination. 

(1) From Example 1, F... = 28,600, 2 dt, = 
0.0807 and = Date fen = 2,311. 

(2) From Fig. 4, for b,,/ts, = 94,¢ = 2.47. From 
Fig. 4, for F.-: = 28,600 Ibs. per sq.in., (1/c)(b./ts,) = 
18.4. Therefore, b, = (1/c)(b./tsn)Ctsn = (18.4)(2.47) X 
(0.025) = 1.135 in., and },/t,, = 1.135/0.025 = 45.5. 
From Fig. 2, for b,/ts, = 45.5, fen = 18,900 Ibs. per 
$q.in. 








(3) 2 Dnln = 0.0807 
(1.135)(0.025) = 0.0284 
New = Dnty = 0.1091 
> Dnata Sn = 2311 
(1.135)(0.025)(18,900) = 536 
New = Dnptn fn = 2,847 


New F,,; = 2,847/0.1091 = 26,100 Ibs. per sq.in. 
(4) From Table 1, for a 24S-T Alclad stiffener 
with vertical legs, K = 1.07. 


F., = KF,.; = (1.07)(26,100) = 27,900 Ibs. per sq.in. 


(5) For this stiffener the corrected test crippling 
stress F,,’ = 28,900 Ibs. per sq.in. (see Table 2), or 
the test stress is 3.7 per cent higher than the predicted 
stress. 


Stiffener Column Curves 


After the crippling stress of a stiffener or stiffener- 
sheet combination has been determined, that stress 
may be used as the cutoff for an approximate column 
curve constructed in the following manner: 


STRESSES 505 

(1) Make a graph with stress as ordinate and L/p as 
abscissa. 

(2) For L/p < 20, make stress equal to F,, of the 
stiffener or combination (except for L/p < 40 for 
24S-T extrusion stiffened panels). 

(3) Construct Euler column curve 


F, = cr*E/(L/p)? 


where 
c = 1.5 for design 
L = length of column, in. 
E = modulus of elasticity, lbs. per sq.in. 
p = radius of gyration of stiffener alone, in. 


(4) Draw a straight line tangent to the Euler curve 
through the crippling stress cutoff point at L/p = 20 
(but at L/p = 40 for 24S-T extrusion stiffened panels). 

Fig. 6 shows a comparison between the predicted 
column curve and the corrected test curve for LS161— 
0.064 24S-RT Alclad stiffeners on 0.040 24S-T Alclad 
sheet. 
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LS 161-0.064 24S-RT Alclad stiffener 
A = 0.150 sq. in. p = 0.481 in. 


Fic. 6, Column curve: 
on 0.040 24S-T Alclad sheet. 


SUMMARY OF TEST RESULTS 


Tables 2, 3, and 4 give a summary of the various 
stiffeners and panels tested as well as a comparison of 
the test and predicted crippling stresses. For L/p ap- 
proximately equal to 20 there were 137 different stiff- 
ener and stiffener-sheet combinations tested. All of 
the specimens were tested in duplicate or triplicate. 
Also tests (in duplicate) were made at other L/p values 
(from one to three additional points) for the panels in 
order to determine column curves. All of the speci- 
mens were tested flat-ended, which accounts for some 
of the scatter. In accordance with general practice, 
the compressive stresses were corrected to minimum 
guaranteed properties using the procedure given in the 
ANC-5 Sec. 1.543; and the test column curves, having 
an end fixity constant c = 3, were corrected to the de- 
sign value of c = 1.5. The tests show the fixity value 
of c = 1.5 particularly to be applicable in the design 
of wing and fuselage compression panels. 

The dimensions of the stiffeners are shown in Table 1. 
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TABLE 3 
24S-RT Alclad Stiffeners, 24S-T Alclad Sheet 
1 2 3 4 5 6 7 8 

Fe’, Lbs. Predicted Margin, 

Stiffener Sheet f per Sq.in., Test F..1, Lbs. per Foe = K Feet, } Pe — 1)100 

Stiffener * Gauge, te: Gauge, ter Corrected Sq.in. K Lbs. per Sq.in (z= ) ; 
LS 160 0.025 25,900 24,100 1.08 26,000 — 0.4 
0.032 0.020 32,500 25,300 1.08 27,300 19.0 
0.032 0.025 28,000 24,400 1.08 26,400 6.1 
0.032 0.032 30,200 24,100 1.08 26,000 16.2 
0.040 35,900 31,800 1.08 34,400 4.4 
0.040 0.020 35,200 29,700 1.08 32,100 9.7 
0.040 0.025 35,100 29,000 1.08 31,400 11.8 
0.040 0.032 30,300 28,000 1.08 30,200 0.3 
0.051 0.020 41,500 33,700 1.08 36,400 14.0 
0.051 0.025 39,800 33,200 1.08 35,800 11.2 
0.051 0.032 38,700 31,500 1.08 33,500 15.5 
0.051 0.040 34,400 30,000 1.08 32,400 6.2 
0.051 0.051 34,000 28,300 1.08 30,600 11.1 
LS 160 Modified 0.032 0.032 37,800 35,500 1.08 38,400 — 1.6 
0.0382 0.032 38,500 31,200 1.08 33,700 14.2 
LS 161 0.025 23,200 22,800 1.08 24,600 — 5.7 
0.032 0.040 23,200 22,800 1.08 24,600 — §.7 
0.040 30,400 29,600 1.08 32,000 — 5.0 
0.040 0.040 28,600 25,400 1.08 27,400 4.4 
0.051 0.025 38,200 32,000 1.08 34,600 10.4 
0.051 0.040 35,600 29,600 1.08 32,000 11.2 
0.051 0.051 32,000 28,000 1.08 30,200 6.0 
0.051 0.064 34,000 26,500 1.08 28,600 18.9 
0.064 49,500 39,600 1.08 42,800 15.7 
0.064 0.040 38,700 36,000 1.08 38,900 — 0.5 
X-LS 165 0.025 25,600 21,700 1.08 23,400 - 9.4 
0.040 34,700 27,600 1.08 29,800 16.4 
0.025 0.025 27,600 27,000 1.08 29,200 — 5.5 
0.040 0.025 33,500 33,500 1.08 36,200 — 7.5 
0.064 0.025 41,000 35,000 1.08 37,800 8.5 
X-785-A 0.025 22,400 19,500 1.08 21,000 6.7 
0.040 31,100 28,000 1.08 30,200 3.0 
0.064 48,500 41,000 1.08 44,300 9.5 











* Sketches shown drawn to scale in Table 1. f No sheet gauge given indicates stiffener tested alone. 


The panels tested had two or three stiffeners, with 
stiffener spacings of 4, 5, or 6 in. in most cases. Many 
of the stiffeners were riveted to the sheet with 1-in. 
spacing, although some were spot-welded. In some 
instances of V-stiffeners (greatest use of spot welds) 
the spots were prestressed as a result of small waves 
left in the stiffener by the forming brake; these waves 
were noticeable only in the heavier gauge sheets. In 
connection with possible evaluation of scatter, it is 
well to note that no definite measurements were made 
of maximum crookedness (initial lateral deflection) to 
length ratio, a factor which the Aluminum Company 
research engineers‘ try to hold to less than 1 part in 
3,000 for stiffeners tested alone. 


DISCUSSION 
It is not the purpose of this discussion to go into any 
long explanations of the reasons for the scatter of indi- 
vidual test points but rather to make a few itemized 
comments concerning (1) the validity of the test results 
as a whole and (2) the use of the summary of results in 


conjunction with the curves for predicting crippling 
stress as an aid in designing stiffened panels. 


Validity of Prediction Curves and Test Results 


Fig. 7 is a statistical presentation of the percentage 
margins of test over predicted stress tabulated in 
Column 8 of Tables 2, 3, and 4. Fig. 7 shows fre- 


quency curves for the various stiffener materials— 
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TABLE 3 (Continued) 
24S-RT Alclad Stiffeners, 24S-T Alclad Sheet 





5 6 7 8 


1 2 3 4 
Fes", LD6. Predicted Margin, 
Stiffener Sheet Tf per Sq.in., Test Fic, Lbs. per Foe = K Fra, Fas" 
Stiffener * Gauge, ts: Gauge, ter Corrected $q.in. K Lbs. per Sq.in (= 7% 1) 100 
X-785-B 0.025 21,600 17,600 1.08 19,000 13:7 
0.040 30,100 25,800 1.08 27,800 8.3 
0.064 41,200 36,700 1.08 39,600 4.0 
X-785-C 0.025 22,000 19,500 0.92 17,900 22.9 
0.040 29,200 29,200 0.92 26,800 9.0 
0.064 43,500 42,200 0.92 38,800 12.1 
X-785-D 0.025 24,300 22,500 1.08 24,300 0 
0.040 36,200 32,100 1.08 34,600 4.6 
0.064 52,700 43,700 1.08 47,200 11.7 
LS 156 0.025 0.025 25,100 22,600 1.08 24,400 2.9 
0.032 0.025 28,200 25,000 1.08 27,000 4.4 
0.032 0.032 27,000 24,000 1.08 25,900 4.2 
0.040 0.025 35,500 28,200 1.08 30,400 16.8 
0.040 0.032 34,400 27,000 1.08 29,200 17.8 
0.040 0.040 33,000 25,800 1.08 27,800 18.7 
LS 157 0.050 0.040 34,500 29,500 1.08 31,800 8.5 
0.051 0.051 30,000 27,800 1.08 30,000 Oo. ° 
0.064 0.032 37,800 33,600 1.08 36,300 4.1 
0.064 0.040 41,500 32,200 1.08 34,800 19.3 
0.064 0.051 39,000 30,600 1.08 33,000 18.2 
LS 152 0.040 39,500 40,900 0.92 37,600 5.1 
0.040 0.040 34,300 35,600 0.92 32,800 4.6 
0.040 0.051 32,900 34,400 0.92 31,600 4.1 
LS 153 0.051 44,100 43,100 0.92 39,600 11.4 
0.051 0.040 39,200 39,400 0.92 36,200 8.3 
0.051 0.051 37,500 37,500 0.92 34,500 8.7 
0.051 0.064 37,000 36,200 0.92 33,300 11.1 
LS 154 0.064 0.051 43,700 40,300 0.92 37,100 17.8 
0.064 0.064 42,600 39,400 0.92 36,200 17.7 
LS 174 0.032 0.025 27,400 27,000 0.92 24,800 10.5 
0.040 0.040 24,100 27,100 0.92 24,900 — 3.2 











* Sketches shown drawn to scale in Table 1. 


curves constructed by indicating on the ordinate the 
number of specimens having margins lying within each 
5 per cent increment on the abscissa and then faired by 
statistical methods to give the normal form. Note that 
there are relatively few extruded stiffener combinations. 
Information from Tables 2, 3, 4, and Fig. 7 is sum- 
marized in Table 5. The positive margins of 6 to 9 
per cent at the peak of the frequency curve indicate 
that if the method of this paper is used for a number of 
stiffeners and panels throughout an airplane, the aver- 
age over strength of the stiffeners and panels will be 
only 6 to 9 per cent. 

Test results show that as much as 6.5 per cent differ- 
ence may occur between the crippling stress of identi- 
cal specimens. On the whole the results appear good; 
however, it is felt the scatter would have been less had 
pin-ended tests been made instead of flat-ended ones 
and had all stiffeners and panels had exactly the same 
L/p for the nominal L/p = 20. Furthermore, it should 
be pointed out that the stiffeners tested alone had no 
lateral restraint and, consequently, failed about the 
axis of least radius of gyration. 


t No sheet gauge given indicates stiffener tested alone. 


Drawing the straight line column curve through F,, 
and L/p = 20 is rather conservative in some cases. 
It is for this reason that making the construction 
through F,, and L/p = 40 is suggested for the extrusion 
stiffened panel curves (note panel); in fact, such a con- 
struction approximates the parabolic form which the 
extrusion stiffened panel test points seem to follow. 


Range of Application as Applied to Stiffener Design 


The solid portions of the curves of Figs. 2 and 3 
represent the range of tests for those curves; the dotted 
portions are extrapolations based on theory. 

A summary of the 6/t and R/t ranges is given in Table 
6. 

The writer feels that the extrapolated portions of the 
curves may be used providing the stiffeners so designed 
have good proportions. The following miscellaneous 
remarks on stiffener design should clarify the preceding 
statement: 

(1) Avoid relatively long outstanding legs. The 
point marked({f) in Table 2 is a good example where 
this rule was not followed, and the stiffener had poor 
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TABLE 4 
24S-T Extruded Shapes, 24S-T Alclad Sheet 
1 2 3 4 5 6 7 8 
F.,’, Lbs. per Predicted Margin, 
Stiffener Sheet f Sq.in., Test Fat Fi. = K Feet, (: ce” Ne 1) 100 
Stiffener* Gauge, tes Gauge, t., Corrected Lbs. per Sq.in. K Lbs. per Sq.in. Free 
Bulb Angles 
LS 222 30,400 27,900 1.08 30,200 0.7 
Double 222 0.020 31,000 25,300 1.08 27,300 13.6 
222 0.020 33,000 26,200 1.08 28,300 16.6 
222 0.025 28,700 24,500 1.08 26,400 8.7 
Double 222 0.025 30,500 25,600 1.08 27,600 10.5 
222 0.032 26,000 23,400 1.08 25,300 2.8 
222 0.040 22,400 22,600 1.08 24,400 — 8.2 
220 31,800 29,500 1.08 31,800 0 
220 0.020 31,700 27,700 1.08 29,900 6.0 
Double 220 0.020 35,000 28,800 1.08 31,100 12.5 
220 0.040 29,100 24,000 1.08 25,900 2.4 
223 31,600 29,500 1.08 31,800 — 0.6 
224 38,000 30,100 1.08 32,500 16.9 
Bulb Tees 
LS 247 29,200 25,500 1.08 27,500 6.2 
248 28,000 25,300 1.08 27,400 2.2 
249 32,200 28,400 1.08 30,600 5.2 
234 40,100 34,000 1.08 36,800 9.0 
Angles 209 30,700 30,600 1.00 30,600 0.3 
200 34,400 33,600 1.00 33,600 2.4 
219 40,100 36,600 1.00 36,600 9.6 
208 32,400 31,900 1.00 31,900 1.6 
210 32,800 31,000 1.00 31,000 5.8 
Channel 264 36,000 32,900 1.00 32,900 9.4 
Tees 236 31,000 28,600 1.00 28,600. 8.4 
232 39,000 35,600 1.00 35,600 9.6 





* Sketches shown drawn to scale in Table 1. 


t No sheet gauge given indicates stiffener tested alone. 


TABLE 5 





Margin at Which —2.5 to +17.5% 
Peak of Normal 


Probable 
Variation 
of Single 
Prediction 


Specimens from 
—2.5 to +12.5% 
Margin, Fig. 7 


Specimens from 


Margin, Fig. 7 














Number* Range* of Test Frequency Curve Num- % Num- % from 
Type Stiffener Combinations Margins, % Occurs, % ber Total ber Total Peak, +% 
24S-T Alclad Sheet 50 —4to +16 5.7 47 94 43 86 3.5 
24S-RT Alclad Sheet 64 —6 to +23 8.5 55 86 43 68 4.6 
24S-T Extrusion 25 —8to +17 6.2 23 92 20 81 4.0 
* Two exceptions of Table 2 excluded. 
TABLE 6 (2) Use vertical-legged stiffeners when possible. 
F Range of Stiffener Tests Table 1 shows them to be more efficient, as would be 
b/t ; —R/t expected. 

Ray ne a he (3) Make stiffener thicker than sheet by one gauge 
fn ae — — ee inet at least when possible. The tests indicate that for pure 
Type Stiffener Free Free Free Free : aacail lof I ff 1 fabri 

24S-T Alclad Sheet 0.9-45.2 3.1-28.9 3.2-10.5 3.6-8.4 C°mpression load a pane of heavy sti leners anc fabric 
24S-RT Alclad Sheet 8.8-42.9 0.6-27.7 3.0-10.6 4.8-8.3 Covering would likely be the most efficient structure on 
24S-T Extrusion 22.0 6.0-19.0 the basis of strength-weight ratio. The point marked 





proportions. A little judgment is the best guide here 
(as elsewhere). This statement is particularly true in 
using the curves of Figs. 2 and 3 with extruded stiff- 
eners of unusual shape and size, since all.the extrusions 
tested had legs of equal thickness with a limited range 
of b/t. 


(**) in Table 2 shows the rule just stated should be 
rigidly followed in the instance of slant-legged stiffeners 
where the angle is considerably off the vertical as is 
that of the LS 174 stiffener. Some of the J-stiffener 
tests indicate that stiffener ‘hinner than sheet by one 
gauge at most may be used for vertical-legged stiffeners, 


although not as efficiently as if the sheet were thinner 














than the stiffener. As stated previously, the tendency of 
the sheet to form waves across the stiffener increases 
with the thickness of the sheet. 


(4) Use flanges on the outstanding legs of stiffeners 
for making the one-edge-free element perpendicular to 
the sheet. This step increases the crippling stress, but 


may be undesirable from the production point of 
view. 


PREDICTING CRIPPLING STRESSES 


CONCLUSION 


If in the future the curves given herein are extended 
and the scatter reduced by making pin-ended tests, it 
may be possible to eliminate many panel tests, which 
should be the case for numerous applications of the 
present method. Even so, it is hoped the article as it 
stands will be of real assistance to structures engineers 
in predicting panel crippling stresses. 
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Rational Analysis of Tension-Field Beams 


H. W. SIBERT* 


Unwersity of Cincinnati 


SUMMARY 

The analysis in this paper is not based upon Wagner’s funda- 
mental assumption that the force in each flange of a tension-field 
beam is in the direction of the flange. The wing is considered as 
a whole, and the stresses in the tension-field shear web are deter- 
mined by a method similar to that for finding the longitudinal 
shear stresses in a beam with a shear-resistant web. Approxi- 
mate formulas are then developed which can be applied to beams 
of either constant or variable cross section. 


INTRODUCTION 


See COMMONLY ACCEPTED THEORY of tension-field 
beams was first given by Wagner! ?* and was 
later summarized by Kuhn.* Wagner’s theory is based 
upon the assumption that the force in each flange of 
a tension-field beam is in the direction of the flange. 
It will be shown later that Wagner’s analysis gives the 
correct solution for a particular type of beam of con- 
stant cross section. For all other beams his analysis 
will give erroneous results. 

Wagner’s fundamental assumption is the result of 
an attempt to simplify the problem of the bending of a 
metal wing by assigning different functions to the 
different parts of the wing—namely, the bending 
moment is to be resisted by the flange material and the 
vertical shear by the shear web in conjunction with the 
vertical components of the forces in the flanges. In this 
paper the wing is considered as a whole and the stresses 
in the tension-field shear web are determined by a 
method similar to that for finding the longitudinal 
beam with a_ shear-resistant 


shear stresses in a 


web. 

This paper will consider only beams with a single 
tension-field shear web. The more complicated prob- 
lem of a beam with more than one shear-carrying 
member will be treated in a subsequent paper. 

WITH A SINGLE TENSION-FIELD 
SHEAR WEB 


CANTILEVER BEAM 


Consider a cantilever beam with a single tension-field 
shear web. Let the x (horizontal) axis pass through 
the center of gravity of the entire flange portion of each 
cross section, and let all the external loads be per- 
pendicular to the x axis. Let A, be the total area of 
the flange portion of a cross section at distance x from 
the free end of the beam. Fig. 1 shows the bending 
stresses f, acting on A, when the beam is bent upward 
by its external loads. The horizontal tensile stresses X 
in the web (due to the diagonal tensile stresses in the 
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direction of the wrinkles) have been replaced by their 
moment M, about the center of gravity of A, and 
their resultant H, acting through this center of 
gravity. 

Let M and M, be the moments about the center of 
gravity of A, of the external loads and of the stresses 
fo in Ay, respectively. M is positive when causing 
compression in the upper fiber, M, is positive when f, 
in the upper fiber is compressive, and M,, is positive 
when the moment of the horizontal tensile stresses X 
on the portion of the web above the center of gravity of 
A, exceeds the moment of X on the portion of the web 
below this c.g. Because of these sign conventions 


M,-M,=M M,=M+M, (1) 


because the two internal moments, M, and My, must 
balance the external moment M. 

In the following analysis all stresses and forces per- 
pendicular to a cross section are positive when tensile. 
Let Ay be the area of the upper portion of A », and let 
Ry be the resultant of the bending stresses f, on Ay. 
Since M, is positive when f, is compressive in the 
upper fiber while f, and Ry are positive when tensile, 
it is evident from Eq. (1) that 


fo = —My/I = —(M + My)y/I (2) 
Sav faA - —(M + My)Qv/I (3) 


where y is the distance of f, from the center of gravity 
of Ay, I is the moment of inertia of Ay, and Qy the 
moment of Ay about this center of gravity. 

Since the resultant tensile force, Hw, on the shear web 
acts through the center of gravity of the flange area 
A», it must be balanced by a direct compressive stress, 
—fa, distributed uniformly over Ap. Let Py be the 


Ry = 
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ANALYSIS OF 
resultant force due to stress fzon Ay. Since fz, Hw and 
Py are positive when tensile, 


fa = —Hy/Ar Py = faAv = —HwAvd/Ar (A) 


Fig. 2 shows the horizontal forces and stresses (in their 
positive directions) acting on a portion of the beam 
bounded by two vertical planes at distances x and x + 
dx from the unsupported end and by a third plane 
through the wrinkle (at angle a with the x axis) which 
passes through the top of the shear web at the left face 
of Fig. 2. As shown in Fig. 2, this wrinkle extends 
downward to the right because the beam was assumed 
to be bent upward by its external loads. Xy is the 
value of X at the intersection of the top of the shear web 
with the left face of Fig. 2. 
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Let ¢ be the thickness of the shear web. Since 
Wagner! has shown that there can be no stress along a 
cut parallel to the wrinkles of a tension-field web, the 
result of equating to zero the sum of the forces in Fig. 2 


in the x direction is 
X ypt(dx tan a + dhy) + dRy + dPy = 0 (5) 


(6) 


X ypt(tan Qa a hy’) — —d(Ry “Le Py) ‘dx 
where iy’ = dh/dx. Replacing Ry and Py by their 


values given in Eqs. (3) and (4) 


tXy(tan a + hy’) = d{[(M + My)Qo/T] + 

HyA U /A F } / dx (7) 
For a beam of constant cross section, iy, Qy/J and 
Ay/A, are all constant and Eq. (7) reduces to 


tXy tan a= [CS + My’)Qv/T) + Hy'A u/Ar (8) 


where S = dM/dx is the total shear force on the cross 
section, My’ = dMy/dx and Hy’ = dH;,/dx. 

By virtue of the sign convention for My and the 
fact that X and Hy are positive when tensile, 


hu hu 
Hy = ; Xtdy My = f yXtdy (9) 
—hL —hL 


where y is the distance of X from the c.g. of Ay, and 
hy and hy are, respectively, the distances from the c.g. 
of A, to the top and bottom of the shear web (Fig. 1). 
Letting X’ = dX /dx, 


TENSION-FIELD 


BEAMS 511 


Hy’ = (10) 


hu hu 
f X'tdy M,’ = [ yX'tdy 

—AL J —AL 
Fig. 3 shows a tension-field shear web of a cantilever 
beam bent upward under the action of external loads. 
Consider a cross section at distance x from the free end 
of the beam. Since all stresses are constant along a 
wrinkle, the horizontal stress XY at any point A on this 
cross section is the same as Xy at the top of the shear 
web (point B) on a second cross section. Thus, the 
determination of Xy at the top of the shear web at each 
cross section by means of Eq. (7) or (8) will define X 
at any point on a desired cross section, except when a 
wrinkle through that point fails to intersect the upper 
boundary of the shear web. However, this situation 
occurs only on cross sections near the free end where the 
web stresses are generally much smaller than at cross 

sections at greater distances from the free end. 
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WAGNER BEAMS 
In this article tension-field beams in which the 


centers of gravity of the upper and lower flanges fall, 
respectively, at the top and bottom of the shear web 
will be called Wagner beams (Wagner! makes this very 
assumption in his analysis of tension-field beams). 
If A, is the area of the lower portion of A, and if the 
moment of inertia of Ay and A, about their own centers 
of gravity are neglected, the following relations hold 
for a Wagner beam © 


hyAy ao hyAz I = hy’?Agy + hy?Arz 
Qu = hyAy (11) 


because fy and hy, are now the distances of the centers 
of gravity of Ay and A;z, respectively, from the center 
of gravity of Ay. By virtue of the relation fyAy = 
h,A, it can be shown that 
I = (hy + hy)hvAu = hQy 

The substitution in Eq. (7) of 1/h for Qy/I and 
h,/h for Ay/Ayr gives an equation from which Xy for 
a tapered Wagner beam can be found for many loading 
conditions. For the sake of simplicity, the actual deter- 


Ap/Ar = h,/h (12) 





512 JOURNAL OF 
mination of Xy for various loading conditions will be 
limited to Wagner beams of constant cross section. 


WAGNER BEAM OF CONSTANT CROSS SECTION 


From Eqs. (8), (10) and (12) the value of Xy for a 
Wagner beam of constant cross section is 


hu 
htXy tana = S+ J 
hb 


Kuhn‘ has found that @ is approximately 45 deg. for 
tension-field beams similar to those found in airplane 
wings. For the sake of simplicity, this value of a will 
be used in the analysis which follows. 


tX’(y + hy,)dy (13) 
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Fig. 4 shows a portion of a Wagner beam of constant 
cross section. The stress X at point A on a cross sec- 
tion at distance x from the free end is equal to Xy at 
point B, whose distance z from the free end is a function 


of y. In fact (Fig. 4), 


y=hy—x+2 dy = dz (14) 


In Eq. (13) let a = 45 deg. and change the variable 


in the integral from y to z by means of Eq. (14). The 
final result is 
htXy = S+ - tXy'(h — x + 2)dz (15) 
x—h 


where Xp’ is the value of dX y/dx with each x replaced 
by a z. Since both x and z are horizontal distances 
from the free end of the beam, this definition of Xy 
makes Xy’dz = dXy. After Xy'dz in Eq. (15) is re- 
placed by dXy, an integration by parts reduces this 
equation to : 


fa tXy ds = § 
x—h 


A change of variable in Eq. (16) from z to y by means 


of Eq. (14) gives 
hu 
f tX dy = S 
—AL 


tT when a = 


(16) 


(17) 


Wagner! shows that X = 45 deg., r 


being the equivalent vertical shear stress in the web. 
Hence, Eq. (17) proves that the entire shear force S on 
a cross section of a Wagner beam of constant cross sec- 
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tion is carried by the shear web. Moreover, the same 
result will be reached when a is not 45 deg. 

Wagner’s assumption that the shear force in each 
flange is parallel to the flange leads to the same con- 
clusion when the beam is of constant cross section. 
Hence, for Wagner beams of constant cross section, 
Wagner’s assumption and the analysis developed in this 
article give identical results. However, the results will 
not be the same for beams other than Wagner beams of 
constant cross section. 

FORMULA FOR BEAMS OF CONSTANT 
Cross SECTION 


APPROXIMATE 


When a beam of constant cross section is not of the 
Wagner type, X must be found from Eq. (8). Since 
Hy’ and My,’ are functions of X’ (Eq. (10)), the deter- 
mination of X from Eq. (8) requires the solution of an 
integral equation. In airplane wings SQ,/J generally 
contributes much more to X than the two remaining 
terms combined. This fortunate situation permits a 
solution of Eq. (8) by successive approximations. 
The resulting formula is comparatively simple if only 
two approximations are used. This formula will be in 
still simpler form if a = 45 deg., which is the approxi- 
mate value of a@ for tension-field beams in airplane 
wings.‘ The first approximation is obtained by letting 
Xy = Xyo in Eq. (8) and putting X’ = 0 in Hy’ and 
My’ (Eq. (10)). Then, with a = 45 deg., Eq. (8) 
becomes , 


tX vo = SQy/I ( l 8) 


From Eq. (18) tXyo’ = wQy/I, where w = dS/dx is the 
external load per unit length of beam (positive when 
upward). 

The second approximation 
tX’ = wQy/I in Py’ and My’ (Eq. (10)). 
45 deg., this changes Eq. (8) into 


tXy = (S+ )Qv/I (19) 


hu hu 
ywdy + (Ay/A rf wdy (20) 
—hAL —AL 


Eqs. (19) and (20) were obtained by replacing tX’ in 
Eq. (10) by wQ,/I. Actually, tXyo’ = wQy/I. How- 
ever, the X’ at any point A on a cross section is the same 
as the X,’ at point B on the upper edge of the shear web 
when AB is a wrinkle through point A (Figs. 3 and 4). 
Hence, the w to use at point A is the load per unit 
length of span at point B (Figs. 3 and 4). Thus, w 
in Eq. (20) varies from the load per unit length of span 
at point C to the load per unit length of span at point 
D (Fig. 3). It will be convenient to call these two 
loads per unit length of span w, and wy, respectively. 
Since the spanwise distance between C and D (Fig. 3) 
is approximately equal to the height of the beam and 
since the ratio of beam height to beam length is small 
in airplane wings, there should be little error in assum 
ing a straight line variation of w between w, and wy. 


results from putting 
With a = 


@ = (Qv/I) 














ANALYSIS OF TENSION-FIELD 


If this is done, Eq. (20) reduces to 


@ = [(wi/2) — (Wmhzt/h)| (h?Qy/T) + 


(wphAy/Ar) (21) 


Wm = (Wy + wz)/2 w, = (2wy + wz)/3 (22 


Thus, w, and w, are, respectively, the loads per unit 
length of span midway between C and D and one-third 
of the way from D to C (Fig. 3). 
APPROXIMATE FORMULA FOR BEAMS OF VARIABLE 
Cross SECTION 


The ratio Ay/Ay in Eq. (7) should be practically 
constant for most airplane wings. If Ay/A, is taken 
constant and a = 45 deg., the application of the for- 
mula for the derivative of the product of two variables 
to Eq. (7) gives 


IXy(1 + hy’) = ((S + Mw’)Qv/I] + 
[((M + My)d(Qp/IT)/dx] + Hw'Av/Ar (23) 


Little error should result from neglecting My in Eq. 
(23) because My is always small numerically compared 
with if the cross section in question is not close to the 
free end of the beam. With My neglected, Eq. (23) 
becomes 


tXy(1 + hy’) = [(S + My’)Qo/I] + 


((Hw’Av/Ar) + Md(Qo/D)/dx] (24) 


The dimensions of Qy and J are cu.in. and in.', 
respectively. Since the height / of an airplane wing 
increases as the distance x from the free end increases, 
Qy/I decreases as x increases and d(Q,/I)/dx is negative. 
Moreover, since M is positive when the free end of the 
beam is bent upward under the action of its external 
loads, Md(Qy/I)/dx will be negative for such a beam. 
Since f increases as x, hy also increases with x and 
hy’ is positive. It would be on the side of safety, then, 
to neglect both Md(Qy/I)/dx and hy’ in Eq. (24). 
If this is done, Eq. (24) becomes 


tXy = [(S + Mw’)Qo/I| + Hw’Av/Ar (25) 


which is identical with Eq. (8). Since Eq. (19) is a 
close approximation of Eq. (8), it is also a close approxi- 
mation of Eq. (25). The neglected terms, hy’ and 
Md(Qy/I)/dx, can now be introduced into Eq. (19) 
in the same places as they were in Eq. (24). Eq. (19) 
then becomes 


tXy(1 + hy’) = [(S + 6)Qu/I] + Md(Qy/I)/dx 


where ¢ is defined by Eq. (21). 

A more accurate approximation of Eq. (24) could 
have been found by letting X’ = 0 in My’ and Hy’ 
(Eq. (10)) to obtain the following first approximation of 
Xv: 


(26) 


tXvo(1 + hy’) = (SQv/T) + Md(Qp/I)/dx (27) 


A second approximation could then be found by putting 
X’ = Xyo' in My’ and Hy’, which would lead to a 
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more accurate equation than Eq. (26). However, 
Eq. (26) is simpler and also conservative, because the 
value tXyo = SQy/I (Eq. (18)) used in the derivation 
of Eqs. (19) and (26) is larger than the value of tX yo 
obtainable from Eq. (27). 

For Wagner beams J/Qy = h (Eq. (12)). This rela- 
tion does not hold in general, but it is probably true 
that for most airplane wings //Q,y is proportional to h, 
or I/Qy = kh, where k is a constant of proportionality. 
By virtue of this relation, 


d(Qy/I)/dx = —h'/kh? = —(h'/h)Qy/I (28) 


where h’ = dh/dx. The substitution of this value of 
d(Qy/I)/dx in Eq. (26) gives 
tXyp(1 + hy’) = (S+ ¢ — Mh’ h)Qv/I (29) 


Since M and h’ are 
Hence, 


where ¢ is defined by Eq. (21). 
positive in airplane wings, Mh’/h is positive. 
the stress Xy at a given cross section of a tapered wing 
is always less than in a beam of the same height whose 
cross section is constant. 

For a Wagner beam J/Qy = h and Ay/Ap = h,/h 
(Eq. (12)). The substitution of these values in Egs. 


(21) and (29) reduces Eq. (29) to 
tXy(1 4. hy’) = (S h) + (wy ‘2) _ (Mh’'/h?) (30) 


For a Wagner beam 
0 and Eq. (30) 


where w; is defined by Eq. (22). 
of constant cross section fy’ = h’ = 
becomes 


tXy = (S/h) + (w/2) (31) 


ACCURACY OF APPROXIMATE FORMULAS 


Before the accuracy of the approximate formulas 
heretofore developed can be checked, the corresponding 
exact solutions must be known. For the sake of sim- 
plicity, only the approximate formula for a Wagner 
beam of constant cross section (Eq. (31)) will be 
checked. The following loading conditions will be con- 
sidered: (a) concentrated load P at end of beam— 
S = P; (b) uniform load—w = k, S = kx; (c) uni- 
formly varying load—w = kx, S = kx*/2; (d) parabolic 
load—w = kx?, S = kx*/3; (e) cubical load—w = 
kx®, S = kx'*/4. 

The following exact values of Xy satisfy both Eq. 
(15) and Eq. (16): (a) concentrated load—tXy = 
P/h; (b) uniform load—tXy = k[(x/h) + ('/2)]; 
(c) uniformly varying—tXy = k[(x?/2h) + (x/2) + 
(h/12)}; (d) parabolic load—tXy = k[(x?/3h) + 
(x?/2) + (hx/6)]; (e) cubical load—tXy = k[(x*/4h) + 
(x3/2) + (hx?/4) — (h'/120)]. 

Since w; in Eq. (31) is the load per unit length of span 
at x — h/3 (Eq. (22)), the following approximate 
values of Xy are obtained from Eq. (31): (a) concen- 
trated load—tXy = P/h; (b) uniform load—tXy = 
(kx/h) + (R/2); (c) uniformly varying—tXy = 
(kx?/2h) + [k(x — h/3)/2]; (d) parabolic load— 
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tXy = (kx?/3h) + [k(x — h/3)?/2]; (e) cubical load— 
tXy = (kx*/4h) + [k(x — h/3)3/2]. 

A comparison of the exact and approximate values of 
Xy shows that Eq. (31) gives no error for a concentrated 
or uniform load, while the error in the other three cases 
are: (c) uniformly varying—kh/4; (d) parabolic load— 
k[(hx/2) — (h?/18)]; (e) cubical load—k[(3hx?/4) — 
(h?x/6) + (11h?/1080)]. 

A comparison of the first term of each error with the 
first term of the corresponding exact solution shows 
that the relative error cannot exceed the following 
values: (c) 0.5(h/x)?; (d) 1.5(h/x)?; (e) 3(h/x)?. 
Hence, toward the root of an airplane wing where h/x 
is small and Xy has its largest values, the relative error 
in using Eq. (31) for these three loadings is extremely 
small. 

The error from using the approximate formulas for 
beams other than Wagner beams could be obtained 
by finding exact solutions from Eqs. (7) and (8). How- 
ever, an estimate of the accuracy of these approximate 
formulas without first finding the exact solution can 
be made in the following manner. The term ¢ in 
Eqs. (19) and (29) is the additional term that is added 
to the first approximation to obtain the final approxi- 
mate solution. The error involved in neglecting all 
successive approximations after the second ought not 


to exceed ¢. Hence, if ¢ turns out to be small com 
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pared with the first approximation, the error itself 
should be small compared with the first approximation. 


TOTAL STRESS IN FLANGES 


With X known at each point on a cross section, Hy 
and M,, for the cross section can be determined from 
Eq. (9). Since the height of an airplane wing is always 
small compared with its length, not much error should 
result in Hy or My if X is assumed to vary as a straight 
line from X, at the bottom to Xy at the top of the shear 


web. With this simplification, Eq. (9) becomes 


Mwy = [(X1/2) — (Xmhz/h)| th? 
(2Xy + X,)/3 


Hy = Xmth 
Xm = (Xup + Xz)/2 XM = 
Note that X,, is the stress at the midpoint of the shear 
web while X,; is the stress one-third down from the 
top of the web. 

After Hy and My have been found from Eq. (32), 
the total stress f, at any point on the flange is (Eqs. (2) 
and (4)) 


(34) 


fe=fht+fhe = —[(M + My)y/I| — Hw/Ar 


There is no need to make any correction for the vertical 
component of the forces in the flanges as is done in Wag- 
ner’s analysis, because all stresses have been included 


in Eq. (34). 
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